Report of the Working Group on Photon and Weak Boson Production 
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This report discusses physics issues which can be addressed in photon and weak boson production in Run II at the Tevatron. 
The current understanding and the potential of Run II to expand our knowledge of direct photon production in hadronic 
collisions is discussed. We explore the prospects for using the IT-boson cross section to measure the integrated luminosity, 
improving the measurement of the W and Z boson transverse momentum distributions, the Z — > bb signal, and the lepton 
angular distribution in W decays. Finally, we consider the prospects for measuring the trilinear gauge boson couplings in 
Run II. 



1. Introduction 

For the next few years, the Fermilab Tevatron Col- 
lider will be the high energy frontier of particle physics. 
The luminosity enhancement provided by the Main In- 
jector will significantly increase the discovery reach of 
the Tevatron experiments over what has been achieved 
with Run I data. It will also move the experimen- 
tal program into a regime of precision hadron collider 
physics. This will make it possible to address open 
questions of high energy physics from several comple- 
mentary directions. 

Understanding the mechanism of production for 
photons, W, and Z-bosons is important for several rea- 
sons. First of all, it provides an opportunity to directly 
test the Standard Model (SM). Second, photon and 
weak boson production often constitute an irreducible 
background in searches for new physics. Finally, a de- 
tailed understanding of the production mechanism for 
these particles is necessary to control the systematic 
errors in precision measurements, such as the determi- 
nation of the W mass (see Ref. [0). 

In this report, we consider several aspects of the 
production of photons, W, and Z bosons which are 
of interest for Run II of the Tevatron Collider. In 
Sec. ||, direct photon production is discussed. A de- 
tailed overview of our current experimental and the- 
oretical understanding of direct photon production in 
hadronic collisions is presented. Direct photon pro- 
duction has long been considered a probe of QCD and 
a source for extracting the gluon distribution of the 
proton. Unfortunately, not all existing (fixed target 
and Tevatron collider) datasets are consistent, and 
our current theoretical understanding of direct pho- 
ton production, especially at small photon transverse 



momenta, is incomplete. Recent theoretical develop- 
ments, however, offer optimism that the long-standing 
difficulties in direct photon production can finally be 
resolved. The enormous number of photon-jet events 
expected in Run II may help to shed light on these 
issues. In particular, photon - jet correlations should 
be helpful in sorting out the source of disagreement 
between theory and experiment. In addition, the kine- 
matic reach in transverse momentum will be greatly 
extended in Run II. 

In Sec. H several important topics associated with 
W and .Z-boson production are discussed. For many 
measurements in Run II, knowledge of the integrated 
luminosity is essential. The integrated luminosity can 
be extracted either from the total inelastic cross section 
or from the cross section of a theoretically well under- 
stood process with high statistics, such as inclusive W 
production. Using the total inelastic cross section to 
determine the integrated luminosity for Run I has led 
to inconsistencies which may well persist in Run II. It 
may thus be advantageous to use the W production 
cross section as an alternative. In Sec. 3.1 we present 



a brief overview of the magnitude of the experimen- 
tal uncertainties in such a measurement. The total 
uncertainty in the W cross section is found to be dom- 
inated by the uncertainty from the parton distribution 
functions, which is considered in more detail in the 
Report of the Working Group on Parton Distribution 
Functions J || . 

In Sec. [$, we also consider the transverse momentum 
distribution of the Z boson, Z — > bb decays, and the 
prospects to measure the lepton angular distribution 
in W decays. The Z px distribution is of interest as a 
test of QCD, and as a tool for reducing uncertainties in 
the transverse momentum distribution of the W . This 
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is important for a precise determination of the mass of 
the W boson (see Rcf . [ Q ) . Searching for a light Higgs 
boson in the range between 110 GeV and 180 GeV is 
one of the prime objectives for Run II. H — > bb decays 
dominate for Higgs boson masses Mjj < 135 GeV. The 
ability to separate the Higgs boson signal in the Wbb 
and ttbb channels from the large QCD background de- 
pends critically on on the bb invariant mass resolution, 
and thus on the measured 6-quark jet energies. Z — ► bb 
decays offer a testing ground for algorithms designed 
to improve the jet energy measurement for b jets and 
are also useful as a calibration tool. The measurement 
of the lepton angular distributions in W decays serves 
as a probe of NLO QCD. The measurement carried 
out by D0 in Run I is statistics limited. While a QCD 
calculation is preferred, large deviations from QCD are 
not excluded. In Run II, this measurement will allow 
for a much improved test of the QCD prediction. 

In Sec. H of this report, we discuss di-boson produc- 
tion. Vector boson pair production provides a sensitive 
ground for direct tests of the trilinear gauge boson 
couplings. A brief overview of the WWV (V = 7, Z), 
ZjV and ZZV couplings is presented and recent ad- 
vances in our theoretical understanding of the NLO 
QCD corrections to di-boson production are described. 
After a brief review of the limits on trilinear couplings 
obtained in Run I, the prospects for strengthening ex- 
isting bounds in Run II are discussed. In addition 
to improving the measurements of WWV and Z"/V 
couplings, it will be possible to determine the ZZV 
couplings via ZZ production with an accuracy of about 
15% in Run II, and to observe the so-called "radiation 
zero" in Wj production. 



2. Direct Photon Production 

The use of direct photon production as an elec- 
tromagnetic probe of hard scattering dynamics has 
a history which covers more than twenty years. As 
in other electromagnetic processes such as lepton pair 
production or deep inelastic scattering, the point-like 
coupling of the photon to charged particles offers some 
simplifications over purely hadronic probes. Compared 
to hadronic jet production, direct photon production 
offers the apparent advantages of having fewer subpro- 
cesses at lowest order and of avoiding the complications 
of jet definitions when measuring or calculating a cross 
section. This latter point means that one can extend 
the range of transverse momenta to smaller values 
for direct photons than for jets. However, in actual 
practice, these apparent simplifications must be tem- 
pered by having to deal with backgrounds from neu- 
tral meson decays, a lower event rate compared to jet 
production, and complications from photons produced 



during jet fragmentation, to name just a few. Never- 
theless, direct photon data provide information which 
complements that obtained from other hard scattering 
processes. Furthermore, photons may be important 
signatures of physics beyond the SM. Therefore, it is 
necessary to understand the "conventional" sources of 
photons before one can fully exploit them in signatures 
designed to look for new physics. 

In this Section, recent work concerning the phe- 
nomenology of initial-state gluon emission in direct- 
photon production in hadron collisions is reported. In 
Sec. 2A, high-mass direct-photon pairs are used to ex- 
plore the impact of such radiation in terms of effective 
parton transverse momenta, fey. At fixed-target ener- 
gies, data on high-p inclusive tt° and 7r°7r° production 
are used to further clarify the arguments presented. 
We then review progress towards fully resummed QCD 
descriptions and present comparisons of a phenomeno- 
logical kr model to recent fixed-target and collider 
data. Possibilities for more extensive studies with 
data from Run II, and the additional information they 
can provide for these considerations, are explored in 
Sec. 2.1.5 . A consistent picture of the observed de- 
viations of NLO perturbative QCD (pQCD) calcula- 
tions from inclusive direct-photon and ir° data is now 
emerging, and we comment on the implications of these 
results for the extraction of the gluon distribution, 
G(x), in Sec. pi. 



Run II has the potential to significantly expand our 
knowledge of direct photon production. Issues related 
to our understanding of the relevant production mech- 
anisms in the kinematic range accessible during Run II 



are reviewed in Sec. 2.2. In this Section, we also discuss 



observables which may help improve our understanding 
of direct photon production as well as experimental 
issue s wh ich can affect the quality of the data. In 
Sec. 2.2.4 some predictions for photon-jet correlations 
are presented. The potential of observables other than 
the usual single photon px distribution to help eluci- 
date the underlying dynamics is discussed. 

Run I data are available on the associated produc- 
tion of a 7 carrying large transverse momentum along 
with a charm quark c whose transverse momentum 
balances a substantial portion of that of the photon [ 
q|. An intriguing possibility is that the data may be 
used to measure the charm quark density in the nu- 
cleon as well as to probe dynamical correlations pre- 
dicted by QCD. These possibilities are discussed in 
Ref. [ ^| where predictions are obtained from a full 
next-to-leading order perturbative QCD calculation of 
p + p -f + c + X at high energy. The associated pro- 
duction of a photon and a heavy quark is not discussed 
further in this report. 
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2.1. Present Status of Direct Photon Produc- 
tion in Hadronic Collisions 1 
2.1.1. Introduction 

Single and double direct-photon production at high 
p have long been viewed as ideal processes for testing 
the formalisms of pQCD. NLO calculations are avail- 
able for both processes [ §, §, f?J §, |). While the 
importance of including gluon emission through the 
resummation formalism was recognized and available 
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for some time for the di-photon process [ 10, |ll| |, it 
is only recently that this approach has been devel- 
oped for inclusive direct-photon production [ |l2[ [H| 
[l4[ |l5| , [l6| , [l7|l . A complete theoretical description 
of the direct-photon process is of special importance 
as it has long been expected to provide one of the 
best measurements of the gluon distribution in the 
proton. The quark-gluon Compton scattering sub- 
process (gq — ► 7(7) shown in Fig. [I] provides a ma- 
jor contribution to inclusive direct-photon production. 
The gluon distribution (G(x)) is relatively well con- 
strained for x < 0.1 by deep-inelastic scattering (DIS) 
and Drell-Yan (DY) data, but less so at larger x [ [ij] . 
Fixed-target direct-photon data can constrain G(x) at 
large x, and consequently has been incorporated in 
several modern global parton distribution analyses [ 
[l9[ |2C| , pi] . More recently, however, both the com- 
pleteness of the theoretical NLO description of the 
direct-photon process, and the consistency of the avail- 
able data sets have been the subject of intense debate [ 
^2], |23|, |24, 25, Z6|, p7[. Direct-photon measurements 
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Figure 1. Leading order diagrams for direct- photon 
production. 

in collider data, and especially the data expected from 
the forthcoming Run II at the Tevatron, provide an 
important testing ground for novel approaches and im- 
provements in the understanding of the direct-photon 
process, and therefore can help resolve the present ar- 
guments. 



1 Contributed by: L. Apanascvich, M. Begel, Y. Ger- 
shtein, J. Huston, S. E. Kuhlmann, D. Partos, J. Womersley, 
M. Zieliriski, and V. Zutshi. 
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Figure 2. Comparison between proton-induced direct- 
photon data and NLO pQCD calculations for several 
experiments as a function of photon xt (= %Pt / \/s) ■ 
(The CDF and D0 data are from Run la, see Figs. p3| 
and p3for the Run lb update.) 



The understanding of single and double direct pho- 
ton yields, and of the more copious high-p ir° pro- 
duction, is of importance for searches for the Higgs 
in the 77 decay mode at the LHC. In addition, Higgs 
production, both at the Tevatron and the LHC, can be 
affected by soft-gluon emission from the initial-state 
partons, and separation of signal and background 
can benefit from a reliable resummation formalism or 
equivalent parton-shower Monte Carlo descriptions [ 

Urns. 

2.1.2. Parton Transverse Momentum 

A pattern of deviation has been observed between 
measured direct-photon cross sections and NLO calcu- 
lations (Fig. ||). The origin of the disagreement has 
been ascribed to the effect of initial-state soft-gluon 
radiation [ [24J] . Correlations between any produced 
high-p particles probe aspects of the hard scatter not 
easily accessible via studies of single inclusive particle 
production. In particular, studies of high-mass pairs of 
particles such as direct photons and tt° 's can be used to 
extract information about the transverse momentum of 
partons, k T , prior to the hard scatter. Whatever the 
source, any transverse momentum between the partons 
will appear as a net p imbalance among the outgoing 
particles produced in the hard scatter, and is there- 
fore reflected in the vector sum of the individual p 
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Vs (GeV) 



Figure 3. (Q ) of pairs of muons, photons, and jets 
produced in hadronic collisions versus yfs. 




z (p T Balance) 

Figure 4. Diphoton Acj), pout, Q T , and z distributions 
for E706 n-Be data at y/s = 31.1 GeV [||. Overlayed 
on the data are the results from NLO [ ^ (dashed) 
and rcsummcd [ [l(J (solid) calculations. PYTHIA [ |7f 
results (dotted) with (k T ) = 1.1 GeV/c are also shown. 



values of the outgoing particles (Q ). If the outgo- 
ing particles are pairs of photons or leptons, then this 
variable should provide a good measure of (k T ) , with 
(k T ) /parton w (Q^) / y2. When the outgoing particles 
are partons, they will hadronize, but the reconstructed 
jets can also yield a measure of (k T ). 

Evidence of significant k T has been found in several 
measurements of dimuon, diphoton, and dijet produc- 
tion; a collection of (Q ) measurements is displayed in 
Fig. | for a wide raageof v^[||§g||§. The 
values of (Q ) are large, and increase with increasing 
y/s. The dijet k T measurements (Fig. ||) agree quali- 
tatively with the dimuon and diphoton results, though 
they have somewhat higher mean values. Such a shift 
is expected since there is also potential for final-state 
soft-gluon emission in dijet events. The values of (k T ) 
per parton indicated by these data are too large to 
be interpreted as due only to the size of the proton. 
From these observations, one can infer that the (k T ) 
per parton is of order 1 GeV/c at fixed-target energies, 
increasing to 3 GeV/c to 4 GeV/c at the Tevatron 
collider, whereas (k T ) would be expected to be of the 
order of 0.3 GeV/c to 0.5 GeV/c based solely on proton 
size. 

The p imbalance between the outgoing particles can 
also be examined using kinematic variables other than 
Q . Given some finite k T , the two outgoing particles 
no longer emerge back-to-back; the azimuthal angle 
between the particles, Acf>, will differ from 180°. The 
transverse momentum normal to the scattering plane, 



Pout, and the p -balance quantity, 

z = ~(f Tl ■ P T2 IV T \) = (P T2 /P T1 ) cos A<t>, (1) 

are two useful variables, each with two possible values 
per pair of objects. 

High-mass direct-photon pairs have been measured 
at the Tevatron [ [|| ||l], [3(| . Distributions as a func- 
tion of A0, pout, Q , and z for such events from 
E706 [ |3l[] are shown in Fig. ^. Overlayed on the data 
are the results from both NLO [ ^ and resummed [ 
|l(| pQCD calculations. There are large differences in 
the predicted shapes. At leading order, each of these 
distributions would consist of a 6 function. While the 
NLO prediction has finite width due to the radiation of 
a single hard gluon, the resummed theory, which also 
includes the effects of multiple soft-gluon emission, is 
in better agreement with the data. This is particularly 
true for Q , where the NLO calculation tends towards 
infinity as Q t — * 0, while the resummed (RESBOS [ |iofl ) 
calculation follows the shape of the data and goes to 
zero. Also shown in Fig. || are the distributions from 
PYTHIA [ |37j, where k T effects are approximated by 
a Gaussian smearing technique. PYTHIA provides a 
reasonable description of the di-photon data using a 
value for (k T ) consistent with the measurements dis- 
played in Fig. ||. Comparisons between CDF and D0 
data lead to similar conclusions [ p9]. There is also 
good agreement between the WA70 di-photon data [ 
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Figure 5. pout distributions for 77, 7r°7r°, 777°, and 
7771- for E706 TrXBe data at = 31.1 GeV Over- 
layed on the diphoton data are the results from NLO [ 
|and resummed [ pQCD calculations. PYTHIA [ 
37) results with (k T ) — 1.1 GeV/c are also shown. 
Overlayed on the tt°tt° and jir are the results from 
LO pQCD calculations [ Q for various values of (fc T ) 
and fixed (q T ) = 0.6 GeV/c. The tt°tt° data have been 
overlayed on the r]ir° data for comparison. 



[}2] and resummed pQCD [[TlJ]. The increased statist- 
ics expected for Run II should allow for more detailed 
comparisons between di-photon data and theory. 
Similar evidence for effects is seen in analyses 



of high-mass 7r°7r°, 7777°, and 771-° pairs by E706 [ |3l|]. 
This is illustrated by Fig. || which shows a comparison 
of the pout distribution from each of these samples. 
The LO pQCD calculation [S, which incorporates k T 
effects using a Gaussian smearing technique similar to 
that used in PYTHIA [ E^jJ, provides a reasonable char- 
acterization of fc T -sensitive variables such as A<j) and 
Pout for (k T ) similar to that measured for di-photons. 
The (k T ) values needed to provide good matches to the 
data for 7r°7r° and 771- are slightly larger than for 77, 
but that is expected since ir° 's emanate from final-state 
quarks and gluons that can produce additional gluon 
radiation. (We use (<7 r )=0.6 GeV/c for the p due 
to fragmentation.) 

2.1.3. k T Phenomenology 

Similar soft-gluon contributions are expected to 
be present in other hard-scattering processes, such 
as inclusive production of jets or direct photons [ 
§§§,§. Resummed pQCD calculations for 
single direct-photon production are currently under 
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Figure 6. Comparison between a threshold resummed 
and a NLO theory calculation for direct-photon pro- 
duction for two scale choices: p /2 and 2p [ p~2| . 
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Figure 7. NLO and NNLO results for direct-photon 
production in hadronic collisions compared to E706 [ 
H pBe data at y/s = 31.6 GeV [ [l|. 



development [ |T|, |D| |§ [TBj, [T§ [[§. Two recent in- 
dependent threshold-resummed pQCD calculations for 
direct photons [ [l^, [l| do not include k T effects, but 
exhibit less dependence on QCD scales than the NLO 
theory (see Figs. || and 0). These threshold-resummed 
calculations agree with the NLO prediction for the 
scale n rs p 12 at low p . and show an enhancement 
in cross section at high p . 

A method for simultaneous treatment of recoil and 
threshold corrections in inclusive single-photon cross 
sections has been developed [ [l7| within the formal- 
ism of collinear factorization. This approach accounts 
explicitly for the recoil from soft radiation in the 
hard-scattering subprocess, and conserves both energy 
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4 6 8 10 12 

Pt (GeV) 

Figure 8. Direct-photon cross section for the E706 
data [ ^5| . The dotted line represents the full NLO cal- 
culation, while the dashed and solid lines respectively 
incorporate pure threshold resummation and joint 
threshold and recoil resummation [ \LTi . 



and transverse momentum for the resummed radiation. 
At moderate p , substantial enhancements from higher- 
order perturbative and power-law non-perturbative 
corrections have been found at fixed-target energies, 
as illustrated in Fig. | in a comparison with the 
E706 direct-photon measurement at ^fs = 31.6 GeV. 
Although the present numerical results are only ex- 
ploratory estimates of the size of expected effects, it is 
already clear that the phenomenological consequences 
are significant. 

Approximate phenomenological approaches to k T - 
smearing have been used in the past to investigate 
deviations between data and NLO pQCD. The un- 
derlying hypothesis is that the leading gluon in NLO 
pQCD calculations is inadequate for describing the 
full initial-state parton shower. Full parton shower 
Monte Carlos such as PYTHIA or HERWIG have been 
used to test this idea at collider energies [ pig . We 
used PYTHIA v5.71 [ ^] and the direct photon process 
to extract the number of initial-state gluons as well as 
amount of net k T present in initial-state gluons, after 
subtracting the gluon with the highest initial state p . 
The number of initial-state gluons is shown in Fig. || 
illustrating that the number is significantly larger than 
the NLO pQCD approximation of either or 1. The 
net k T present, after subtracting the highest p gluon, 
is shown in Fig. for 10 GeV/c and 50 GeV/c direct 
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Figure 9. Number of initial-state gluons in PYTHIA 
from the direct-photon process at y/s — 1.8 TeV. 
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Figure 10. Net p due to all but the leading gluon in 
PYTHIA for direct-photon events. 



photons in pp collisions at y/s = 1.8 TeV. The net p 
of such remnant gluons is 2.6 GeV/c for direct photons 
with p = 10 GeV/c, and 5.2 GeV/c for p t = 50 GeV/c. 

At fixed-target energies, parton-showering models do 
not provide sufficient smearing because shower devel- 
opment is constrained by cut-off parameters that en- 
sure the perturbative nature of the process. Since tra- 
ditional NLO calculations do not account for the effects 
of multiple soft-gluon emission, a kinematical model 
was employed to incorporate k T effects in available 
pQCD calculations of direct-photon (and tt°) produc- 
tion [ ^J. The relationship between this phenomeno- 
logical fc T -smearing and the Collins-Soper-Sterman 
(CSS) resummation formalism [ Q was considered 
in some detail in Ref . [ E!J . 



6 




"t3 10 J - 

a 



p T (GeV/c) 

Figure 11. The variation of k T enhancements, K(p ), 
for the E706 pBe data at V? = 31.6 GeV. 



The same LO pQCD [ |3f| program that successfully 
characterized high-mass pair production was used to 
generate K-factors (ratios of LO calculations for any 
given (k T ) to the result for fc T =0) for inclusive cross 
sections (Fig. |ll]). These K-factors were then applied 
to the NLO calculations. This procedure involves a risk 
of double-counting since some of the /c T -enhancements 
may already be contained in the NLO calculation. 
However, the effects of such double-counting are ex- 
pected to be small [ Q . 

As illustrated in the upper part of Fig. [n], the 
K-factors for direct-photon production at E706 are 
large over the full range of p , and have p -dependent 
shapes — a behavior reminiscent of that obtained from 
the full resummation formalism [ [TtJ . The lower part of 
Fig. [ll] displays K-factors for 7r° production, based on 
the same model. The data appear to require somewhat 
larger values of (k T ) in the case of 7r°s. At the same 
values of (k T ) the K-factors in 7r° production are some- 
what smaller than for 7 production. This is expected 
because 7r°'s originate from the fragmentation of the 
partons. 

Figure [ll] also displays the K- factor for photons used 
by the MRST group [ |2J in recent fits to parton distri- 
butions. Their result was obtained through a different 
technique involving analytical smearing of the param- 
eterized photon cross section, rather than an explicit 
parton-level calculation. Although the correction is of 
similar size, it has a different p -dependence. It should 
be noted, however, that despite the similarity in the 
values of the K-factors used by the MRST group and 
the ones presented here, the (k T ) values cited by MRST 
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Figure 12. The photon and ir° cross sections for E706 
pBe data at ^/s = 31.6 GeV [ Q compared to the 
fc^-enhanced NLO calculations. 



are lower by a factor of 2. This difference can be traced 
to an erroneous relation between parameters of their 
analytical smearing function and the transverse mo- 
menta of partons in the hard scatter [|2l[| . This has al- 
ready been pointed out in Refs. [|24|, |27[. The approach 
used here [ Q is based on explicit parton kinematics 
and therefore does not suffer from this problem. 

The treatment of fc T -enhancements proposed in 
Ref. [ , based on parton distributions unintegrated 
over the parton transverse momenta, suggests possible 
modifications of the above simple picture. Reference [ 
|27| imposes strong ordering of momentum transfers of 
emitted gluons, which prevents transverse momenta of 
the incoming partonic system from approaching p ; k T 
values are correlated with the scale at which the par- 
ton distributions are sampled. In their approach, the 
K-factors are expected to be smaller than those shown 
in Fig. H|, and have less p -dependence. Additional 
scrutiny of the theoretical ideas should help resolve 
these differences. 

2.1.4. High-p Production 
Fixed Target 

Invariant cross sections for inclusive direct-photon 
and 7T° production are displayed for the E706 pBe 
data at y^£= 31.6 GeV in Fig. |lj, with overlays from 
theory [ [45[ . Discrepancies between NLO pQCD the- 
ory (dotted curves) and the data are striking. The 
enhancements, generated using (k T ) values consistent 
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Figure 13. Comparison between the E706 direct- 
photon data and NLO pQCD calculations with and 
without k T enhancements, for several values of (k T ). 



Figure 15. Comparison between proton- induced ir° 
data and NLO pQCD calculations for several exper- 
iments as a function of xt- 
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Figure 14. Comparison between the WA70 and UA6 
direct-photon data and NLO pQCD calculations with 
and without k T enhancements, for several values of 
(k T ). 



with the data on high-mass pairs (Fig. ||), can ac- 
commodate both the shapes and normalizations of 
direct-photon and ir° inclusive cross sections. 

Comparisons between direct-photon data from 
E706 [ ||, WA70 [ H, pi, and UA6 [ @ and 
/c T -enhanced NLO pQCD [ [24] are shown in Figs. [l3] 
and 14. The values of (k T ) were based on the data for 
high-mass pairs from E706 (Fig. ||), and the di-photon 
data from WA70 [f§|,f| ((k T ) = 0.9±0.1±0.2 GeV/c). 
The center-of-mass energies for WA70 and UA6 (y/s rs 
24 GeV) are smaller than those for E706. Correspond- 
ingly, (k T ) values for these experiments are expected to 
be slightly smaller than the values required for E706 
(Fig-!). 

A recent survey of 7r° production found that cur- 
rent NLO pQCD calculations significantly undershoot 
the data [ |2(|. A comparison between 7r° data and 
NLO pQCD is shown in Fig. [l5| for several experi- 
ments 
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The data are consistently a factor of 2 
to 3 above theory. The above phenomcnological model 
should also be valid for pion production. Using (k T ) 
similar to, but slightly higher than that for direct 
photons, good agreement is obtained for 7r°'s mea- 
sured by E706 [ @, WA70 [ and UA6 [ || 
(Figs. [l6| and [l?]). The fc T -enhanced predictions com- 
pare well with the ir° cross sections, with all the E706 
and UA6 direct-photon data, and with the 7r _ beam 
direct-photon cross sections of WA70. 

Since fc T -smearing affects similarly direct-photon 
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Figure 16. Comparison between the E706 ir° data and 
NLO pQCD calculations with and without k T enhance- 
ments, for several values of (k T ). 
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Figure 18. Ratios of direct-photon cross sections to 

the 7T cross sections as a function of p for various 

_ T 

experiments at several values of yfs. 



• <k T > = 0.0 GeV/c 
■ <k T > = 0.9 GeV/c 
o (k T >= 1.1 GeV/c 



7t"p -> Jl"X (WA70) 

Vs=23.0 GeV 
- -0.45 < x„ < 0.55 



....... 1 












pp->lt°X(UA6) 






Vs=24.3 GeV 




-0.1 < <0.9 





pp n"X (WA70) 

Vs=23.0 GeV 
-0.35 < x„ < 0.45 





O • ? « * 



pp -> ji X (UA6) 
Vs=24.3 GeV 
-0.1 <y<0.9 



P T (GcV/c) 



Figure 17. Comparison between the WA70 and UA6 7r° 
data and an NLO pQCD calculation with and without 
k T enhancements, for several values of (k T ). 



and 7T° data, the ratio of direct-photon to ir° pro- 
duction should be relatively insensitive to k T . Ex- 
perimental and theoretical uncertainties also tend to 
cancel in such ratios. Figure [l8] shows the ratios of 
cross-sections [ |5(J for direct-photon to ir° production 
for WA70, UA6, E706, and R806 [ || The 
results from WA70 (Fig. |l|) and UA6 (Fig~g), at 
approximately the same y/s, appear to differ signifi- 
cantly The ratio for NLO theory differs from that of 
WA70 by a factor of three, but only by » 30% from 
the UA6 data. The WA70 and UA6 tt° results agree 
(Fig. |l|) and most of the difference is therefore in the 
direct-photon cross section. Similar j/ir a comparisons 
are shown for E706 at y/s = 38.8 GeV (Fig. and 
R806 at y/s = 63 GeV (Fig. The same 30% level 
of agreement can also be found for the E706 data at 
•v/s = 31.6 GeV and the R806 data at yfi = 31 GeV 
and 45 GeV. 

The discussion of direct-photon data from fixed- 
target and ISR experiments [ while rejecting the 
k T interpretation for the observed deviations, pointed 
to limitations for the applicability of NLO calculations 
at lower-p values and to inconsistencies among experi- 
ments used in comparisons to NLO theory. Given the 
shortcomings of a purely NLO description, a global 
comparison of the available data to fully resummed 
theory may provide useful insights into these issues. 
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Figure 19. 7/71"° comparison for the WA70 pp data at 
^/s = 23 GeV. Overlayed are the results from the NLO 
pQCD calculations. 



Figure 21. 7/71"° comparison for the E706 pBe data 
at y/s = 38.8 GeV. Overlayed are the results from the 
NLO pQCD calculations. 




Figure 20. 7/71" comparison for the UA6 pp data at 
y/s = 24.3 GeV. Overlayed are the results from the 
NLO pQCD calculations. 



Figure 22. 7/71° comparison for the R806 pp data at 
y^i = 63 GeV. Overlayed are the results from the NLO 
pQCD calculations. 
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Figure 23. A comparison of the CDF isolated direct- 
photon data at y/s — 1.8 TeV from Run lb with a 
NLO pQCD prediction and two implementations of 
soft gluon corrections to the NLO prediction [ Q . 
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Figure 25. A comparison between CDF and UA2 data 
at yfs = 630 GeV. 
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Figure 24. A comparison of the D0 isolated direct- 
photon data at y/s — 1.8 TeV from Run lb with a 
NLO pQCD calculation in both the central and for- 
ward rapidity regions [ p7| ■ 



Run I Results 

The consequences of k T smearing are expected to 
depend on yfs (Fig. ||). At the Tevatron collider, where 
p is large compared to k T , the above model of soft- 
gluon radiation leads to a relatively small modification 
of the NLO cross section. Only the lowest end of the 
P t spectrum is modified significantly, and the K-factor 
exhibits the expected ~ 1/p 2 behavior for a power 
correction. 

A comparison of the results from Run lb, shown 
in Figs. [2^ and [24| [ |56|, |57j, confirms the expected 
deviation in shape at low p . Using di-photons, CDF 
has measured (k T ) = 3.6 ± 0.8 GeV/c at \/s = 
1.8 TeV [ |33j. Employing this value, the phenomeno- 
logical model adequately describes the shape of the 
data in Fig. ^3|. Also shown in this figure is an im- 
plementation of soft-gluon corrections using an en- 
hanced parton shower 
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The phenomenological-fc T 
model provides a better agreement with data than is 
available in the enhanced parton shower model. The 
agreement between the phenomenological model imple- 
mentation of k T smearing and the direct-photon data 
can also be seen in preliminary results from CDF at 
y/a = 630 GeV (Fig. ||). 

The CDF data in Fig. |2^ have been normalized up- 
wards by a factor of 1.2 for the benefit of a shape com- 
parison. Without this normalization, then the CDF 
data lie below the NLO pQCD calculation at high p . 
Figure also contains a comparison between CDF 
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Figure 26. Left: Distributions of azimuthal angle dif- 
ference A(j) between the direct photon and a recoil 
charged track for several values of k T and for two 
photon p ranges. Right: Similar dependences for the 
photon pout- 



and UA2 [ g |§ data at y/s = 630 GeV, where a 
similar deficit is observed at high p for the UA2 data. 
There is currently no explanation for this effect. 

2.1.5. Expectations for Run II 
fc„ Studies at Low-n 

To study the level of k T induced by multiple gluon 
emissions at the collider, the experiments have em- 
ployed the relatively low-statistics di-photon data. In 
Run II, both CDF and D0 will have precision mag- 
netic tracking, which will permit studies of k T effects 
using "two-arm" data on pairs consisting of a direct 
photon and a recoil charged track, in a spirit similar to 
that of the fixed-target investigations of jit and 7r 7r° 
pairs discussed before. The advantage of this approach 
is that it will obviate the need for jet reconstruction 
(difficult at low p ), and minimize complications from 
jet energy scale calibration. Of course, photon+jet 
systems are also of interest, but may be harder to study 
in the range of interest for checking the effects of k T 
smearing. 

We simulated the expected behavior of photon-track 
systems using the same LO Monte Carlo employed 
in previous fixed-target studies. The results are illus- 
trated in Fig. ^6| for Acj) and pout of the photon, for 
two representative ranges of photon p that span the 
region where the k T effects appear to Be important in 



the inclusive photon cross sections from Run I. (For the 
measurement oipouT for photon, the scattering plane 
is defined by the colliding beams and the recoil track.) 
The sensitivity to the value of (k T ) in the range of 
to 3.5 GeV/c is clearly seen for both variables. The A(j> 
distributions become narrower with increasing photon 
p^, as expected from simple kinematic arguments. For 
a fixed (k T ), the width of the pout distributions is 
relatively insensitive to p , and can therefore be par- 
ticularly useful for mapping out the dependence of k T 
on event kinematics (especially on the p of the pho- 
ton). To properly interpret the widths of such distri- 
butions in terms of k T induced by gluon-radiation, it 
is important to subtract the amount generated in the 
fragmentation of partons into the charged particles (as 
was done in the fixed-target analyses). 

Additional handles on interpreting the data can 
be obtained through studies of p distributions of 
charged particles from fragmentation of partons recoil- 
ing against direct-photon triggers. In the presence of 
significant initial-state k T , such distributions are ex- 
pected to become softer than expected from standard 
fragmentation functions (determined e.g. from e + e~ 
data), since the fc T -kick tends to increase the photon 
p while taking it away from the recoil side. This should 
be observed most clearly at low photon p 's, where the 
effect from k T is greatest. 

Photon Purity 

The particularly interesting connection between two- 
arm studies and measurements of the inclusive photon 
cross section is mainly at relatively low p (10 GeV/c 
to 35 GeV/c), where Run I results indicate significant 
deviation of the cross section from expectations from 
NLO pQCD (Figs. || and ||). This is also the region of 
high statistics, and, consequently, where detailed stud- 
ies will be possible. It is therefore important to achieve 
a high photon-signal purity in this region, which has 
proved to be difficult in Run I. 

In the case of D0, calorimeter response will be mod- 
ified in Run II by the presence of a central solenoid 
magnet and preshower detectors [ ^|. The separa- 
tion of photon signal and background in Run I was 
based on the fraction of electromagnetic energy de- 
tected in the first longitudinal layer of the calorimeter 
(-Eemi/^tot)- This quantity is particularly sensitive 
to differences in the early stages of shower develop- 
ment initiated by single photons, and by photons from 
decays of 7r°'s of the same p . In the Run II configu- 
ration, the signal-background discrimination based on 
this variable is expected to deteriorate, but simulations 
indicate [ ^HJ that Run I performance can be recovered 
by using instead the fraction of energy deposited in 
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Figure 27. Comparisons of simulated discriminations 
between photons (solid) and 7r°'s (open) for the D0 
detector in the Run I configuration (top, discriminant 
variable is log(i?EMi/£ , TOT); for the same discriminant 
in the Run II configuration (middle); and for the Run II 
configuration and log^cps/^TOT) discriminant (bot- 
tom). 



Figure 28. Additional 7 vs ir° discrimination with the 
D0 detector in Run II is possible using differences in 
the respective distributions for the number of recon- 
structed preshower clusters matched with the calorime- 
ter shower (left) or for the width of the preshower 
cluster (right). Solid lines are for 7's, dashed for 7r°'s. 



the preshower detectors (-Ecps/^tot), as illustrated 
in Fig. [27] for photons at central rapidities. Using 
the fine-grained shower-profile information from the 
preshower strips (ca. 7 mm triangles), an additional 
factor of two rejection of background (while maintain- 
ing high signal-efficiencies) has been achieved in our 
Monte Carlo simulations [ |6l[] . As illustrated in Fig. |28|, 
multiple preshower clusters for sufficiently large sep- 
aration between photons from meson decays can be 
resolved, or inferred from the widths of clusters when 
the showers are not fully separated. Results are shown 
for p « 15 GeV/c and 7/ ss 0.95, where the separation 
of photon signal from background is particularly diffi- 
cult. Another variable useful for the discrimination is 
the distance in the R-(j> plane between the detected 
CPS cluster and the photon position (the latter cal- 
culated from the primary vertex and the position of 
the calorimeter cluster) [ |S2[ . Figure ^ shows this 
distribution for p w 15 GeV/c and r\ w 0.1. 

Even higher background rejection can be achieved 
by exploiting correlations between the calorimeter and 
preshower shower profiles in a covariance-matrix ap- 
proach. A study was undertaken using a simplified 
version of the Run I covariance matrix, with added 
preshower variables [ p2[. Figure p(j shows the Et 



dependence of the 7r° efficiency for % 2 cutoffs corre- 
sponding to 90% and 60% photon efficiency. 

In the central region, CDF tools for 7-background 
separation in Run I (shower width in the electro- 
magnetic shower-maximum detector, and conversion in 
the central-preradiator detector) will remain the same. 
These tools provide a clear separation of the photon 
signal and the 7r°-dominated background. The addi- 
tion of a new scintillator-based endplug calorimeter 
with a preshower and shower-max detector will offer 
an extension of these tools to the forward region in 
Run II. 

Thus, we expect a better signal purity at low p in 
Run II than was achieved in Run I, which will facilitate 
more precise measurements of low-p direct photons 
and di-photons. 

Photons at High-p 

At high values of p , separating the direct-photon 
signal from background and the expected reach in p 
are determined by the collected luminosity. As shown 
in Sec. 2.2, for the initial luminosity of 2 fb -1 , the 
inclusive direct-photon cross section measurement can 
be extended beyond 300 GeV/c. While the high-p 
data will permit detailed tests of perturbative QCD, it 
is not expected to be sensitive to differences in recent 
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Figure 29. The residual of CPS cluster position relative 
to the photon or 7r° in the plane in the D0 detec- 
tor. Clusters from 7r°'s have larger residuals, especially 
when one of the decay photons does not convert in the 
solenoid. 
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Figure 30. Efficiency for a tt° to pass covariance matrix 
X 2 cuts corresponding to photon efficiencies of 90% 
and 60%. (Preliminary study for the D0 detector in 
Run II.) 



paramcterizations of gluon distributions at large x (eg., 
CTEQHJ and CTEQM). Thus, the determination of 
the large- a; gluon distribution will have to continue to 
rely on direct-photon data from fixed target experi- 
ments and on jet data at the Tevatron. 

2.1.6. Impact on the Gluon Distribution 

The largest uncertainty in any parton distribution 
function (PDF) is that for the gluon distribution. At 
low x, the gluon can be determined indirectly from 
scaling violations in quark distributions, but a di- 
rect measurement is required at moderate to large 
x. Direct-photon production has long been regarded 
as potentially the most useful source of information 
on the large-x part of the gluon distribution. And 
direct-photon data, especially from CERN fixed tar- 
get experiment WA70, have been used in several 
global analyses [ ^l], |63| . Another process sensitive to 
the gluon distribution, through the gluon-gluon and 
gluon-quark scattering subprocesses, is jet production 
in hadron-hadron collisions. Precise data from Run I 
are available over a wide range of transverse energy 
and, indeed, provide a constraint on the gluon distri- 
bution in an x range from about 0.05 to 0.25. How- 
ever, the low statistical power of the jet cross section 
at high Et, and the dominance of the qq scattering 
subprocess in that kinematic region, do not provide 
for a similar constraint at large x. 

Figure |3l| sh ows several gluon parton distribution 
functions [ |18| that provide a reasonable fit to the 
DIS and DY data used in the CTEQ4 fits [ ||. 
The excursions shown (normalized to the CTEQ4M 



gluon) provide an estimate of the uncertainty in the 
gluon distribution. The gluon distribution seems 
reasonably well-constrained by these data, except at 
large x. Also shown in Fi g. |3l| are the gluon dis- 
tributions from CTEQ5HJ [ |64| (which best fits the 
Tevatron high-Er jet data) and three recent MRST 
PDF's [ HI. The MRST PDF's incorporate the 
WA70 direct-photon data using the fc r -enhancements 
described above. Nevertheless, the theoretical prob- 
lems associated with direct-photon production have 
discouraged the CTEQ collaboration from using the 
direct-photon data in their recent fits [ Q . The recent 
work on resummation offers hope that this situation 
will change. 

The CTEQ4HJ PDF's were determined by increas- 
ing the weight for the CDF jet cross section at 
high Et [ jl^ ]. In the resultant fit, the increase was 
achieved through a significant increase in the gluon 
contribution at large x, without inducing serious con- 
flicts with any of the other experiments used in the 
CTEQ4 data sets. This increase was allowed by the 
uncertainty in the gluon distribution in this x range, 
a flexibility not present for any of the quark distri- 
butions. Another demonstration of the uncertainty 
in the gluon distribution at large x can be seen in 
Fig. H|, where the CTEQ4HJ gluon distribution is 
plotted along with that of CTEQ4M and three recent 
MRST PDF's. At x fa 0.6, there is over an order 
of magnitude spread between the CTEQ4HJ and the 
MRST g I gluon distributions. 

The CTEQ4HJ gluon distribution (and its successor, 
CTEQ5HJ) provides the best description of not only 
the CDF jet cross section, but that of D0 as well. At 
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Figure 31. The ratio of gluon distributions obtained 
in the study in Ref. [ [if] to the CTEQ4M gluon at 
two different Q scales. All of these gluon distributions 
correspond to PDF's which provide a reasonably good 
fit to the CTEQ4 DIS/DY data set. 



x 0.5, corresponding roughly to Et « 450 GeV for 
central r], a doubling of the gluon distribution (com- 
pared to CTEQ4M) results in only a 20% increase in 
the inclusive jet cross section. On the other hand, 
the fixed-target direct-photon yield of E706, produced 
mainly through gq scattering, is proportional to the 
gluon distribution. A fit using the E706 direct-photon 
data and the k T K-factors is also shown in Fig. |32|; 
this result is very similar to the CTEQ4M gluon. At 
the highest reach of the E706 data, the CTEQ4HJ 
gluon distribution is a factor of 4 to 5 larger than the 
CTEQ4M gluon. With the advent of more complete 
theoretical treatments 
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tion, the E706 data should have great impact on the 
determination of the behavior of the gluon at large x. 
This would have implications not only for fixed-target 
direct-photon data, but also for collider physics at 
highest Q 2 scales. 

2.1.7. Conclusions 

Direct-photon physics remains a viable and interest- 
ing program for Run II of the Tevatron. The Run II 
measurements of single and double direct photons, and 
of photon and jet or single track correlations, will reach 
larger p and have improved detection efficiency at low 
p , compared to Run I. Although the data are not ex- 
pected to improve directly our knowledge of the gluon 
distribution at intermediate and large x, it can do 
so by providing a testing ground for newly developed 
theoretical models and formalisms, and by helping 
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Figure 32. A comparison of the CTEQ4M, MRST, and 
CTEQ4HJ gluons, and the gluon distribution derived 
from fits that use E706 data [ H}. The g | and g j. 
gluon densities correspond to the maximum variation 
in (k T ) that MRST allowed in their fits. 



clarify the currently confused role that multiple gluon 
emission play in direct-photon production (and other 
high-p processes). Once this physics is properly un- 
derstood, the existing fixed-target data should provide 
one of the best constraints on the gluon distribution, 
as has been envisioned for a long time. 

We have examined the best available experimental 
information on production of single and double direct 
photons (and mesons) at large p in both fixed-target 
and collider energy regimes. Recent theoretical de- 
velopments offer optimism that the long-standing dif- 
ficulties in the proper description of these processes 
can finally be resolved. While there is still no final 
consensus, the trend of recent developments has led 
to an increased appreciation of the importance of the 
effect of multiple gluon emission, and to the emergence 
of tools for clarifying this issue. 

To summarize, measurements of the production of 
high-mass pairs of high-p particles at WA70, E706, 
CDF, and D0 provide consistent evidence for the 
presence of large k T . NLO pQCD calculations [ ||, 
which include effects due to the radiation of a sin- 
gle hard gluon, compare poorly to fc T -sensitive dis- 
tributions in di-photon data. RESB0S [ 0, a NLO 
pQCD calculation, which also includes the effects of 
multiple soft-gluon emission through the CSS resum- 
mation technique, compares well with the shape of 
the di-photon data. LO pQCD calculations [ B7L B8| 
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that incorporate k T effects through Gaussian smear- 
ing techniques, provide reasonable characterizations of 
distributions for pairs of direct photons and mesons. 

While the apparent inconsistencies between differ- 
ent direct-photon and 7r° data sets are not under- 
stood [ ^4], [2^, ^6) , we found it instructive to consider 
results on the 7/71- ratio from WA70, UA6, E706, and 
R806. Various experimental and theoretical uncer- 
tainties tend to cancel in such a ratio, which is also 
relatively insensitive to fc T -effects. We find that the 
ratio from theory agrees to w 30% with data from 
UA6, E706, and R806 over the range 24 GeV < ^ < 
63 GeV. 

LO pQCD has been used to estimate the impact of 
k T on the inclusive production of high-p^ direct photons 
and 7r°'s. This simple phenomenological model is able 
to account for differences between NLO pQCD calcula- 
tions and inclusive data over a wide range in ^/s. While 
the approximate nature of such models is clear, and has 
been discussed in several recent papers, the emerging 
formalism for the full (threshold and recoil) resumma- 
tion of inclusive direct-photon cross sections appears 
to vindicate much of the understanding of effects from 
multiple gluon emission that has been achieved using 
approximate tools. The resummation formalism can be 
expected to provide a solid foundation for the treat- 
ment of k T , at which time a global reexamination of 
parton distributions, with an emphasis on the determi- 
nation of the gluon distribution from the direct-photon 
data, should become possible [ ^35| . 

2.2. Direct Photon Production at the 

Tevatron 1 
2.2.1. Introduction 

In this section we summarize features of direct pho- 
ton production which are relevant in the kinematic 
range to be covered in Run II. In Sec. 2.2.2] , a compari- 
son between Run I data and the corresponding theoret- 
ical description is presented. Several potential problem 
areas are noted. Sec. 2.2.3 contains a brief description 



of the kinematic reach expected for Run II, based on 
an integrated luminosity of 2 fb _1 . Included here is 
a discussion of the sensitivity to parton distribution 
functions and to what extent direct photons at Run II 



can help constrain the gluon distribution. In Sec. 2.2.4 



some predictions for photon-jet correlations are pre- 
sented. The potential of observables other than the 
usual single photon pt distribution to help elucidate 
the underlying dynamics is also discussed. 

2.2.2. Comparison to Run I Data 

Data for the inclusive cross section for direct photon 
are available over a wide range in energies from fixed 
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target and collider experiments. By now it is well 
known that it has not yet been possible to simulta- 
neously describe all of the experimental results with a 
next-to-leading-order (NLO) QCD calculation. A pat- 
tern of discrepancies between theory and experiment 
exists in both the fixed target and the collider data sets. 
This situation has been reviewed in [ 



. An 

analysis similar to that for direct photon production 
in Ref. [ |25| has been performed for the case of 7r° 
production [ [2(| . The two processes are closely related 
since 7r°'s decaying to two photons provide much of 
the background which must be dealt with when ex- 
tracting the signal for direct photon produ ctio n. Some 
of these issues are also dealt with in Sec. 2.1 where a 



detailed comparison to fixed target data is also pre- 
sented. While it is clearly of interest to understand 
direct photon production over the entire range of avail- 
able energies, this section will focus on those aspects of 
the data which can most directly be addressed during 
Run II. The first step is to examine the theoretical 
description of the data from Run I. 

A comparison of NLO QCD predictions to the di- 
rect photon data from CDF and D0 has indicated 
the presence of a deviation of the data from theory 
at low values of transverse momentum [ ^3], [24]]. This 
deviation decreases if the effects of soft gluon radia- 
tion are taken into account by applying a Gaussian 
kx smearing model using a value of (kr) measured in 
di-photon production in the two experiments.^] Such a 
kx treatment is phenomenologically motivated. 

Recently, there has been progress in more sophisti- 
cated treatments of soft gluon radiation near threshold 
in the parton-parton scattering process [ [l4| [36], 
At large values of transverse momentum for the pho- 
ton, the phase space for the emission of additional glu- 
ons in the hard scattering is limited. This limitation 
on the emission of real gluons upsets the balance in 
the theoretical expressions between virtual and real 
emission contributions. The result is large logarithmic 
corrections near the threshold for the parton-parton 
scattering subprocesses. These large corrections can 
be resummed in a relatively compact formalism. The 
results [ |l2|, 13 indicate that the corrections to 
existing next-to-leading-order calculations are large as 
xt — ^Pr/y/s approaches 1. Away from the re- 
gion at the edge of phase space it is observed that 
the corrections to the NLO results coming from the 
threshold resummation are relatively small over much 
of the xt range covered in the fixed target experiments 
when the renormalization and factorization scales are 
chosen to be pt/2 and the resummed results show 
an overall reduction in the sensitivity to the choices 



1 Contributed by: J. Huston and J.F. Owens 



*The value of kx can be directly measured in di-photon events 
since the photon 4-vectors can be measured precisely. 
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of these scales. However, the threshold resummation 
corrections alone are not sufficient to explain the dis- 
crepancies observed between the theoretical predic- 
tions and some of the fixed target experimental results. 
In addition, threshold resummation cannot explain the 
deviations observed by D0 and CDF at the low px end 
of the measured distributions. 

Another approach to resumming soft gluon effects is 
that of Ref. [ |7| which uses the DDT [ |]J or g T -space 
method. This technique has recently been applied to 
vector boson production [ |6^] and compared with the 
impact parameter method of Ref. [ Q . In Ref. [ |2?j 
a parton-parton luminosity function is defined which 
depends on the net transverse momentum of the pair of 
colliding partons. The parton distributions are probed 
not at the scale of the hard scattering process, but 
instead at a scale given by qx- This means that the 
scale of the parton distributions is typically somewhat 
smaller than, for example, the transverse momentum of 
the produced photon. This results in an enhancement 
of the cross section. In addition, the parton transverse 
momentum is taken into account when the final photon 
Px is calculated. Some enhancement of the fixed target 
predictions is noted using this technique, but it is in- 
sufficient to fully describe all of the fixed target data. 
For the collider energy range, Ref. [ quotes only 
a small effect, which does not appear to be sufficient 
to explain all the observed deviations. At present, the 
resummation is done only with the leading-log terms 
included. More results from this technique are antici- 
pated as next-to-leading-log terms are included as well. 

Quite recently a new formalism for simultaneously 
incorporating both the threshold resummation and the 
resummation of kx (or recoil) effects has been devel- 
oped [ 0. This formalism possesses the desirable 
property of simultaneously conserving both energy and 
momentum in the resummation process. The initial 
results presented in [ |l7j indicate that the threshold 
enhancement referred to above is correctly reproduced 
while, in addition, there is a large enhancement due 
to the newly included recoil effects. Detailed appli- 
cations and studies of the scale dependence and the 
dependence on non-perturbative input parameters are 
expected soon. 

The initial indications of discrepancies between the- 
ory and experiment in the collider data came from 
results taken in Run la. A comparison of the higher 
statistics results from CDF in Run lb, shown in Fig. |23|, 
confirms the shape deviation at low px and the agree- 
ment with the kx smearing correction obtained using 
the Gaussian smearing model. (Also shown is an at- 
tempt to implement the soft gluon corrections using an 
enhanced parton shower [ |22j. In this case, the effect 
seems to drop off more slowly than the deviation ob- 



served in the data.) Note that the kx smearing correc- 
tion falls off roughly as l/p T , as expected from a power 
correction type of effect. This is not true in the case 



of fixed target experiments, as discussed in Sec. 2.1 



where the steeply falling parton distribution functions 
enhance the effects of the soft gluon radiation. 

The D0 direct photon cross section for the central 
region is also consistent with such a deviation at low 



Pt [ 57 , as shown in Fig. |24| while no conclusion can be 
reached for the D0 cross section in the forward region 
(see the lower half of Fig. |J and Ref. [ |7|). The 
possibility that the discrepancy between theory and 
experiment may be dependent on rapidity is interesting 
and is one that can be investigated in more detail with 
the higher statistics expected from Run II. 

In an attempt to achieve a better theoretical descrip- 
tion of the data it is important to investigate what flex- 
ibility exists within the conventional QCD hard scat- 
tering NLO formalism. In order to perform perturba- 
tive QCD calculations, one must specify the renormal- 
ization and factorization scales. For the latter, there 
arc two scales, corresponding to the factorization of 
collinear singularities for the initial state parton dis- 
tributions and the final state fragmentation functions. 
For most hard scattering calculations these three scales 
are chosen to be proportional to the characteristic large 
transverse momentum - that of the photon in this case. 
Often, all three scales are set equal to each other. How- 
ever, this latter step is not necessary and it is reason- 
able to ask whether or not one can describe the shapes 
of the CDF and D0 data in the region below px = 30 
GeV by a suitable variation of all three of the scales [ 
p9[ . In Fig. [53| the Run lb CDF data are compared to 
several curves corresponding to different choices of the 
renormalization and factorization scales. One can see 
that it is possible to get a steepening of the slope at the 
low-pT end, but only at the price of an increase in the 
overall normalization. Apparently, it is not possible to 
get both the shape and the normalization correct by 
such a strategy. 

In Fig. |23|, the CDF data has been normalized up- 
wards by a factor of 1.2 for an easier shape comparison. 
If this normalization is taken out, as shown in Fig. [54|, 
it is evident that the data fall below the NLO QCD 
prediction at high px- Also shown is a comparison to 
the CDF photon data taken at y/s — 630 GeV and the 
photon data from UA2 [ [5^] where a similar deficit is 
observed at high px- For most observables, typically, 
the data lie above the NLO QCD predictions so this is 
somewhat of an unusual situation. It is interesting to 
note that a measurement of the photon fraction (of the 
photon candidate sample) indicates that the photon 
fraction seems to be leveling off at approximately 80%, 
rather than saturating the sample at near 100% at high 
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Figure 33. Comparison of the CDF Run lb data with 
several NLO curves corresponding to different choices 
of the renormalization and factorization scales. 



Figure 34. A comparison of the CDF photon data from 
Run lb (at both ^ = 630 GeV and ^ = 1800 GeV) 
with a NLO QCD prediction and the direct photon 
data from UA2. 



Pt- The latter outcome would be predicted from very 
basic considerations: with a fixed Et isolation cut, 
jets are required to fragment into a ir° with a higher 
z value as the pt of the photon candidate increases. 
Such a fragmentation is suppressed by the sharp falloff 
of the fragmentation function at high z. It will be 
interesting and important to understand this behavior. 
The increased statistics of Run II will allow both the 
low pt and high pt regions to be investigated more 
thoroughly. 

2.2.3. Expectations for Run II 

For the purposes of this section we shall assume an 
integrated luminosity for Run II of 2 fb _1 . As noted in 



Sec. 2.2.2, the data on direct photon production from 



Run I extend to a transverse momentum of approx- 
imately 120 GeV. The increased statistics expected 
from Run II greatly extend this range as shown in 
Figs. |35| and |3(| These figures have been generated 
using the next-to-leading-logarithm program of [ |7^] 
with the CTEQ5M [ |i| parton distributions and with 
the renormalization and factorization scales set equal 
to pt/2. The errors shown are statistical only and 
the results are presented for the transverse momen- 
tum range where more than 10 events are expected 
in a 10 GeV bin of pt with a total integrated lu- 
minosity of 2 fb _1 and 4 fb _1 , respectively. No effi- 
ciency/acceptance corrections, however, have been ap- 
plied to these estimates. Typically, these corrections 



are on the order of 50%. 

One can see that the range of useful statistics 
extends past 300 GeV in the first case and past 
350 GeV in the second. This extended coverage in pt 
corresponds to an increased range of sensitivity for the 
values of the parton momentum fractions of the collid- 
ing hadrons. Recall that the inclusive cross section for 
single photon production in the central region is sensi- 
tive to average values of the parton momentum fraction 
x approximately equal to xt = ^Pt/^/s. A range of 
xt out to about 0.3 will be covered, corresponding 
to a similar range of (x) for the parton distributions. 
This extended range suggests that the relative ratios 
of the underlying subprocesses should change signifi- 
cantly over the pt range to be covered. To investigate 
this, the inclusive cross section is displayed in Fig. [[j] 
along with the contributions from the various parton 
scattering subprocesses. For ease of comparison, the 
same results are shown in Fig. |3^ on a linear scale 
relative to the total rate. The results in both of these 
figures were generated using the leading-logarithm ap- 
proximation with the CTEQ5L parton distributions. 
The leading-log approximation was chosen so that the 
separation between the point-like and fragmentation 
contributions would be unambiguous. Processes with 
two or more partons in the final state in addition 
to the photon can populate regions of phase space 
where the photon is collinear with one of the par- 
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Figure 35. Direct photon px distribution for Run II 
with errors based on an integrated luminosity of 2 fb -1 . 
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Figure 37. Contributions of the various subpro- 
cesses for direct photon production calculated in the 
leading-log approximation using the CTEQ5L parton 
distributions. 
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Figure 36. Direct photon px distribution for Run II 
with errors based on an integrated luminosity of 4 fb -1 . 



tons. These topologies correspond to the fragmenta- 
tion process. Thus, the higher order terms mix the 
contributions from the point-like and fragmentation 
(or bremsstrahlung) components. One must define, 
through the use of appropriate experimental cuts, pre- 
cisely what is meant by the bremsstrahlung compo- 
nent. For the purposes of the discussion being pre- 
sented here the leading- log predictions are sufficient. 
Note that the magnitude of the next-to-leading-order 
corrections relative to the leading-log predictions for 
the choice of scale used here (pr/2) is slowly varying 
over the kinematic region being studied, as shown in 
Fig. H 

The results shown in Figs. |7| and [58] illustrate sev- 
eral points worth noting. First, the fragmentation 
component is expected to be a negligible fraction of 
the total rate. This is due primarily to the imposition 
of an isolation cut which rejects events with more than 
1 GeV of hadronic energy accompanying the photon in 
a cone of radius 0.4 about the photon direction^] Such 
a cut is necessary experimentally in order to control the 
copious background to the photon signal from jets frag- 
menting into 7r°'s. This cone isolation energy is almost 
completely saturated by the underlying event energy 
accompanying the hard scatter, leaving little room for 



^This cut was used for the CDF Run lb direct photon mea- 
surements. The value of the isolation energy cut in Run II will 
have to be increased somewhat due to the contributions to the 
isolation energy from the larger number of minimum bias events 
expected in the same crossing. 
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Figure 38. The same results as shown in Fig. |37] except 
on a linear scale and normalized to the total rate at 
each value of pr- 



energy from the fragmentation of the jet. Note that the 
precise value of the fragmentation contribution relative 
to the total rate will vary as one includes higher order 
effects, but the overall contribution is still expected to 
be small. Next, one sees the dominance of the qg — > jq 
subprocess in the region out to about 100 GeV in pt, 
i.e., the range covered by the Run I data. Beyond 
this range the qq — > jg subprocess becomes dominant. 
Gluon-gluon initiated processes are not expected to 
play a significant role over the px range shown. 

One of the classic applications of direct photon pro- 
duction is to provide constraints on the gluon distribu- 
tion in global fits of parton distributions. The gluon 
distribution is especially uncertain in the region be- 
yond x ~ 0.15 [ [l8|]. Run I results on high-p^ jet 
production from the CDF [ fQ) and D0 [ Q collab- 
orations favor a gluon distribution which is larger at 
high- a; than was anticipated from global fits which did 
not emphasize the high-pr jet data. One such example 
is the CTEQ5HJ [ £4[ set of distributions which are 
favored by both sets of jet data. One might hope that 
the direct photon data could shed some light on this 
issue, but such is not expected to be the case. As shown 
in Fig. the ratio of the CTEQ5HJ and CTEQ5M 
predictions is consistent with unity within about 5% 
over the pt range under consideration. This can be 
understood by referring back to Fig. |3^ where it is 
shown that the gq subprocess decreases in importance 
precisely where one would like to gain a constraint on 



Figure 39. NLO/LO ratio for direct photon produc- 
tion. 



the gluon distribution. 

2.2.4. New Measurements 

As noted in the introduction, the process of direct 
photon production is of interest for a variety of rea- 
sons. In particular, it offers a probe of hard scattering 
dynamics which complements that of jet production. 
At this point, the comparison between theory and ex- 
periment has not reached the quantitative level that 
one would like, so there is clearly more to be done. 
Furthermore, direct photon production has the poten- 
tial to provide constraints on the parton distribution 
functions, particularly that of the gluon. Nevertheless, 
the less than satisfactory description of the current 
data has meant that this role is, as yet, unfulfilled. 
Finally, one must understand conventional sources of 
energetic photons before being able to confidently use 
photons as part of the signal for new physics. With 
these points in mind, it is natural to ask whether there 
are additional observables which could help shed light 
on some of the theoretical problems. 

To date, almost all of the experimental results for 
direct photon production have been for the transverse 
momentum distribution of the photon, i.e., the single 
photon inclusive cross section. Clearly, there is addi- 
tional information to be gained by studying the joint 
rapidity distributions of the photon and a jet, e.g., 
da / dpTdy^dyjet where pt represents the transverse 
momentum of the photon. This observable has been 
measured by the CDF Collaboration [ (73) for photons 
in the the central region with 16 GeV < pr < 40 GeV. 
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Figure 40. Ratio of the direct photon cross sections 
calculated with CTEQ5HJ and CTEQ5M with errors 
based on an integrated luminosity of 2 fb _1 . 
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Figure 42. Jet rapidity distribution for photon in a 
forward rapidity region calculated with the CTEQ5L, 
CTEQ5M, and CTEQ5HJ parton distributions with 
the program of Ref . [ [TOl . 
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Figure 41. Jet rapidity distribution for photons in the 
central rapidity region calculated with the CTEQ5L, 
CTEQ5M, and CTEQ5HJ parton distributions with 
the program of Ref. [ [7(j . 



In this case the theoretical description of the jet rapid- 
ity distribution is good. However, no information is 
yet available for other values of photon rapidity. Some 
leading-log and next-to-leading-log predictions for this 
distribution are shown in Fig. ^ for central photons 
and Fig. for forward photons. Notice the character- 
istic broadening of the jet rapidity distribution as the 
photon is moved forward. Note, too, the change in the 
ratio between the leading-log and next-to-leading-log 
predictions as the photon rapidity is varied. Predic- 
tions are shown for both the CTEQ5M and CTEQ5HJ 
parton distributions. There is only a slight difference 
between the two sets, due primarily to the fact that the 
curves were generated for px > 10 GeV. The resulting 
low values of Xt result in the parton distributions being 
probed in a region where there is little difference be- 
tween the CTEQ5M and CTEQ5HJ sets. On the other 
hand, that means that for these distributions there 
will be relatively less uncertainty due to the parton 
distributions and that, therefore, such measurements 
may be useful in helping to understand the px re- 
gion where there is a discrepancy between the exist- 
ing collider data and the theoretical results. In this 
regard the D0 results for the inclusive photon yield 
shown in Fig. ^ are interesting as they may indicate 
that the theory/data discrepancy has some rapidity 
dependence. The lesson here is that photon-jet joint 
observables would be helpful in sorting out the source 
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of the disagreement between theory and experiment. 
Note that if one wanted to increase the sensitivity to 
differences between parton distributions, then a larger 
minimum pt cut could be employed. 

there is some indication 



As discussed in Sec. 2.2.2 



that the photon yield above a transverse momentum of 
about 30 GeV may actually be less than the theoretical 
predictions, in contrast to the situation at lower values 
of pt- It was mentioned that this may be related to the 
observed behavior of the j/tt ratio. This issue is com- 
plicated by the necessity of placing isolation cuts on 
the electromagnetic triggers in order to reduce the 7r° 
background. NLO Monte Carlo programs can simulate 
the effect of these isolation cuts, but the experimental 
dependence on the parameters of the cuts has not yet 
been compared to that of the theory. If the theoreti- 
cal treatment of the isolation cuts is wrong, then the 
comparison of the NLO results to collider data must be 
considered suspect. In recent years the question of the 
theoretical treatment of such isolation cuts has been 
studied by several authors [ [74| |75|, [7(J . What is needed 
is a data set showing how the cross section depends on 
the cone size and the energy threshold utilized in the 
isolation cuts. Similar comparisons for jet production 
have helped refine the various jet algorithms and have 
been very useful in understanding issues related to the 
theoretical description of jets. 

2.2.5. Conclusions 

Run II offers many opportunities to refine our under- 
standing of the production of photons in hard scatter- 
ing processes. The kinematic reach in transverse mo- 
mentum should be greatly extended and the statistical 
precision of the data will also be increased. While there 
are still problems with the theoretical description of the 
existing data, the Run II data have the potential to 
shed some light on these issues. In particular, data for 
photon-jet correlations and for the dependence of the 
cross section on the parameters of the photon isolation 
cuts will be helpful. A better theoretical understanding 
of direct photon production will enable this process to 
be better used in the study of large transverse momen- 
tum processes and the search for new physics. 

3. Topics in Weak Boson Production^ 

The very large number of W and Z boson events 
CDF and D0 will collect will yield precision mea- 
surements of the W mass and width [ [l], which are 
fundamental parameters of the Standard Model and 
thus need to be determined with the highest possible 
accuracy. The Run II vector boson datasets will of 
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course provide other important advances in the field 
of electroweak physics, and will be the starting point 
of most new physics searches; but they will also be- 
come standard tools for the understanding of many 
sources of systematic uncertainty in otherwise unre- 
lated physics studies. In this Section we will give an 
overview of some of the additional studies it will be 
possible to carry out with vector bosons in Run II; for 
rare decays of W bosons and other topics not covered 
here see Rcf . [ [fTj . 

3.1. The W Cross Section as a Luminosity 
Monitor 

At both D0 and CDF, the uncertainty on the total 
integrated luminosity is approximately 4.5% [ [78| fr9| |. 
Though measuring the luminosity is difficult to do well, 
any measurement at the collider that has an absolute 
normalization depends on it. In several important 
cross section measurements made during Run I, this 
uncertainty contributed greatly to the overall measure- 
ment uncertainty. In particular, for both of the W 
and Z boson cases, the uncertainty in the luminosity 
measurement far outweighed the other systematic un- 
certainties in the measurement (see Table |l|) . 

There is also a continuing controversy regarding the 
actual value of the luminosity at the two detectors at 
the Tevatron [ [7§ At both D0 and CDF, the 

integrated luminosity is normalized by the total inelas- 
tic cross section in pp collisions. CDF has made this 
measurement [ Q and uses it to normalize their lu- 
minosity. D0 did not measure the total inelastic cross 
section, and chooses to normalize to the world average, 
which is ~ 6.2% higher than the CDF measurement [ 
|7Sj |. This is primarily due to a 2.8 standard deviation 
disagreement between the CDF measurement and the 
E811/E710 measurements [ |l[ §2[ g ||. 

In the context of Run II, there are some general 
concerns regarding the measurement of the total inte- 
grated luminosity, particularly at high instantaneous 
luminosity. New luminosity monitors are being in- 
stalled at each detectors as part of the modifications for 
Run II and CDF plans to repeat the measurement of 
the total inelastic cross section. However, it is not clear 
that this will result in a more precise determination of 
jC, and disagreement regarding the measured pp cross 
section is likely to persist. Given the high precision 
of the current W and Z cross section measurements, 
increasing precision in the matrix clement calculations 
(now at NNLO [ ^U), and the expected abundance of 
W and Z bosons produced during Run II (~ 10 6 W 
bosons and ~ 10 5 Z bosons within the fiducial volumes 
of each detector), we may consider using the rate of W 
boson production to normalize the integrated luminos- 
ity. In this subsection, we present a brief overview 
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Table 1 



Measurements of a(W) at the Tevatron. The three 
quoted uncertainties are statistical, systematic, and 
luminosity. 



Detector 


Channel 


<t{W) ■ B (nb) 


D0(1A) 
D0(1A) 
D0(1B) 
D0(1B) 
CDF(IA) 


W ^ ev 
W->fjt,v 
W -> ev 
W ->tv 
W ^ ev 


2.28 ±0.02 ±0.08 ±0.10 
2.02 ±0.06 ±0.22 ±0.09 
2.31 ±0.01 ±0.05 ±0.10 
2.22 ±0.09 ±0.10 ±0.10 
2.49 ±0.02 ±0.08 ±0.09 



of the magnitude of the experimental uncertainties in 
the integrated luminosity using W boson production, 
and discuss the requirements for such a measurement 
to be competitive with standard luminosity measure- 
ments and in what context it would surpass the current 
precision. 

3.1.1. a(W)/C in Run I/II at D0 

As an example, we describe some details of how the 
Run I measurement of the cross section for W boson 
production at D0 would translate into a determination 
of the total integrated luminosity. 

Reversing the relationship between the integrated lu- 
minosity (£) and the production cross section (a(W)) 
in the calculation done by D0 [ J78|, we have the fol- 
lowing relation: 

~ _ N cand (l — Jqcd) — Nz ,„.. 
eA(l + ^)a(W) U 

where N canc i is the number of W boson candidates 
observed, Jqcd is the fraction of candidates expected 
from QCD multi-jet production, Nz is the number of 
candidates that are Z bosons in which one of the elec- 
trons was unobserved, e is the event identification effi- 
ciency, A is the geometric and fiducial acceptance, A T 
is the acceptance times branching ratio for W — » tv, 
and a(W) is the predicted cross section times branch- 
ing ratio for producing W bosons. 

Table |^ show the values of each of the quantities as 
measured in Run I, and the fractional uncertainty each 
would contribute to a measurement of C. The total 
uncertainty from measured quantities alone is 2.2%. 
Including a 4% uncertainty in the prediction of a(W) 
increases the total uncertainty to 4.6%. 

We note two things: 1) even the Run I measurement 
results in an uncertainty on C that is competitive with 
the directly measured resultjj] and 2) by far, the domi- 
nant uncertainty is in the prediction of a(W). 

+Of course, we cannot reliably use the theory prediction without 
the experimental confirmation from Run I! 



Table 2 

Uncertainties on the components of the measurement 
of C if one used the components of the Run I mea- 
surement of the W production cross section by D0. 



Component 


Value 


Error on C 


Iqcd 


0.064±0.014 


1.5% 


e 


0.671±0.009 


1.3% 


A 


0.465±0.004 


0.9% 


-*V cand 


67078 


0.4% 


N z 


621±155 


0.3% 


A T / A 


0.0211±0.0021 


0.2% 


a(W) 


22.2±0.9 nb 


4% 



In preparation for Run II, both D0 and CDF are 
undergoing major upgrades to the detectors. For D0, 
this includes the addition of a solenoid magnet and 
a complete replacement of the tracking system. The 
total luminosity is expected to increase by a factor of 
~ 10. We now consider how each of the factors in Eq. || 
will be affected in Run II. 

The background fraction for events from QCD multi- 
jet production is dominated by systematics for the de- 
tection of electrons. We do not expect the uncertainty 
to decrease dramatically, nor do we expect the overall 
background level to change much. It may be possible 
to reduce the backgrounds and make them easier to 
understand if we use the muon channel; however, no 
serious study has been made on the subject. 

The uncertainty on the lepton identification effi- 
ciency is equal parts Z boson statistics and background 
subtraction statistics and will scale by ~ 1/VlO (we 
assume an integrated cross section of 1 fb _1 in this 
section) . We expect the tracking efficiency to increase, 
but the efficiency of the isolation and $! T require- 
ments may decrease due to multiple interactions and 
decreased resolution due to the presence of the solenoid 
in the Run II detector. Speculating, these effect may 
balance each other, leaving an overall reduction in the 
uncertainty by 1/3 to 0.9%. 

The geometrical and fiducial acceptance will stay 
approximately the same. The uncertainty is dominated 
by the electromagnetic energy scale (0.00319 of 0.004) 
whose uncertainty will be smaller in Run II due to an 
increased number of Z boson events and extra han- 
dles provided by a central magnetic field in the D0 
detector. Scaling the rest of the uncertainties with the 
luminosity, we should be able to halve the uncertainty 
on the acceptance to 0.5%. 

The uncertainties in the expected number of Z can- 
didates and the acceptance of electrons and muons 
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Table 3 

Estimated values for the uncertainties in a measure- 
ment of C in Run II. 



Component 


Error on C 


Iqcd 


1.5% 


e 


0.9% 


A 


0.5% 


cand 


0.08% 


N z 


0.0% 


A T /A 


0.0% 


a(W) 


4% 



from W — ► tv are dominated by MC statistics and 
can be shrunk to a negligible value. 

Finally, the number of candidates will increase by 
a factor of 10 just from the increased luminosity, by 
approximately a factor 1.12 from the increase of the 
center of mass energy from 1.8 TeV to 2.0 TeV, and 
by another factor of ~ 2 if one includes the muon 
channel, resulting in about 1.4 x 10 6 candidates, which 
translates into a statistical uncertainty of about 0.08%. 

Table | summarizes the expected experimental un- 
certainties for measuring C in Run II at D0. A mea- 
surement of the luminosity will be dominated by the 
background level and uncertainty and by the uncer- 
tainty in the prediction of the cross section. Figure [l3| 
shows how changes in the QCD background level and 
uncertainty, and uncertainty in the cross section pre- 
diction, affect the overall uncertainty in the resulting 
luminosity. In each case, all quantities from experi- 
ment were fixed to their Run I values except the one 
being varied. The uncertainty on the cross section 
was made negligible when considering the sensitivity 
to the background level and uncertainty. We note 
that decreasing the fractional uncertainty on the QCD 
multi-jet background by a factor of 2 - 4 reduces the 
luminosity uncertainty from 2.2% to 1.7 - 1.8% - a 
factor of 0.3 — 0.2. The relatively large uncertainty in 
the background level (20%) translates into a significant 
sensitivity to the background level itself. Essentially, 
the trade-off for loosening the selection criteria and 
allowing more background into the data sample is that 
the background must be measured much more precisely 
in order to maintain a small uncertainty. In Run II 
(as in Run I), the best strategy to minimize the ef- 
fect of the background uncertainty on the total cross 
section measurement will be to minimize the level of 
the background altogether - if the background can be 
made negligible, the effect of the uncertainty of the 
background on the cross section (no matter how large) 
will also be negligible. The cross section uncertainty 



continues to dominate the situation. However, if the 
cross section uncertainty is kept to approximately the 
size of the experimental uncertainty, the resulting lu- 
minosity uncertainty is more than a factor of two better 
than the Run I value. 
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Figure 43. Impact of changing the overall background 
level (b), the uncertainty on the background level (Sb), 
and the predicted a(W) in determining C. Except 
when being varied, the experimental values are kept 
to the Run I determinations. Except when varying the 
uncertainty on the cross section, Sa(W) is set to 0. The 
circles denote f = 0.05,0.10,0.15,0.25. The squares 
denote background fractions, b = 5, 10, 15%. The 
triangles denote = 0.001,0.005,0.01,0.03,0.05. 

The solid line shows the current values of the luminos- 
ity uncertainty at D0 and the dashed line shows the 
uncertainty if the rate of W boson production was used 
and a(W) was known perfectly. 



3.1.2. Counting on the Prediction for a(W) 

If one determines the integrated luminosity using the 
rate for W boson production, there are two fundamen- 
tal issues that must be resolved. First, one must decide 
that the calculation for a(W) is reliable in itself; that it 
agrees with experiment. Second, one must determine 
the uncertainty in the calculation, since it will likely 
dominate the uncertainty in C. 

The first issue is likely not problematic. If all cross 
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section measurements were normalized to a specific 
a(W) calculation, then the worst impact of a change 
in that calculation would be to modify all measure- 
ments in the same manner. Additionally, the advan- 
tage would be that all measurements would be easy to 
compare, since the controversy over the measurement 
of the total inelastic cross section of pp collisions would 
be circumvented. 

The problem of determining the uncertainty on the 
calculation is more difficult. Various ad hoc methods 
have estimated the uncertainty on a(W) to be 3 — 5% [ 
[78).The uncertainty in the cross section is dominated 
(almost exclusively) by uncertainties in the PDF's 
which go into the calculation. (The uncertainties due 
to higher order QCD and electroweak corrections are 
likely much smaller than these.) Two efforts to un- 
derstand the uncertainties in the PDF's quantitatively 
are described in the report of the Working Group on 
Parton Distribution Functions [ || . 

3.1.3. Conclusion 

We have described the current expectations for mea- 
suring the integrated luminosity in Run II using the 
rate for W boson production. The experimental un- 
certainties, totaling ~ 2.2%, are dominated by the 
uncertainty in the background from QCD multi-jet 
interactions. However, the total uncertainty is dom- 
inated by an ill-defined uncertainty on the prediction 
for the production cross section for W bosons. With 
sufficient progress in the continuing effort to quantify 
this uncertainty, we may be able to reliably determine 
the total integrated luminosity in Run II using the rate 
for W boson production to ~ 3 — 4%. 

3.2. Determination of the Weak Boson pt Pro- 
duction Spectrum 

During Run II at the Tevatron, CDF and D0 will 
obtain the largest data set of e+e - pairs resulting from 
via pp — ► 7*/Z to date, pushing the analysis of vector 
boson production characteristics over the edge from be- 
ing limited by statistical uncertainties, to being limited 
by systematic uncertainties. The di-electron final state 
provides two important experimental handles. Elec- 
trons themselves are among the best-measured objects 
at either detector, with far better resolution in en- 
ergy and position than most final-state high-pT objects 
(such as jets or muons). Additionally, the di-electron 
final state provides complete kinematic information 
about the hard collision; the four-momentum of the 
Z/'j* state is known unambiguously. Since the elec- 
troweak character of the decay to di-electrons is gener- 
ally uncorrelated with the QCD characteristics of the 
production of the vector boson state, di-electron pro- 
duction in pp collisions via the Drell-Yan production 
is a sensitive probe for investigating many aspects of 



QCD. The rapidity (y) distribution is the Drell-Yan 
analog to deep inelastic scattering structure function, 
providing additional constraints on PDF's. The trans- 
verse momentum (pr) distribution is sensitive to pre- 
dictions from standard perturbative QCD at high-px 
(~ Q 2 ) and to predictions from soft-gluon resumma- 
tion calculations at low-px [ ■ Additionally, the 
low-px region is sensitive to non-perturbative effects 
not calculable in pQCD. These effects are included 
via a universal form factor, not unlike PDF's, whose 
parameters must tuned to data. 

In addition to the intrinsic benefits of precision mea- 
surements of QCD, there are practical benefits for 
other measurements at the Tevatron. In the low-px 
region, where the cross section is highest, uncertainties 
in the phenomenology of vector boson production have 
contributed to the uncertainty in the measurement of 
the mass of the W boson (Mw) Diboson, top 

quark, and Higgs boson production all have single and 
di-electron backgrounds from W and Z boson produc- 
tion that will be more constrained through a precise 
measurement of Z/j* production properties. Also, the 
universality of the resummation approach requires fur- 
ther experimental testing, with implications ranging 
from the impact on the precise determination of Mw, 
to the production of Higgs bosons and di-photons [ 

High mass Drell-Yan e + e pairs are experimentally 
distinctive: the electrons typically have large fir's, are 
separated from each other in rj and cf>, and tend to 
be separated from jets and other activity in an ev- 
ent. In Run I, CDF and D0 collected such events 
with electrons in the central (|f?dct| < 1.1) and forward 
(1.1 < |r?det| < 2.4/2.5) regions, providing a coverage 
in the e + e~-pair rapidity of up to \y\ ~ 3. Since 
the online and off-line electron identification efficien- 
cies are comparable (for both experiments), additional 
off-line requirements primarily enhance the rejection 
of background from QCD multijet events that were 
mis-measured as electrons. At CDF, the most power- 
ful discrimination is provided by tracking and precise 
electron shower centroid measurements. To cope with 
the large jet backgrounds in the forward detector re- 
gion, CDF developed the SVX-Plug tracker [ and 
improved the matching of tracks in the vertex tracker 
with calorimeter shower positions during Run I. Apply- 
ing such forward tracking techniques to both e's of a 
pair reduced the backgrounds from about 10% to a per- 
cent or less. At D0, the most background rejection was 
obtained by additional isolation and shower-shape re- 
quirements, reducing the backgrounds from ~ 10—15% 
to 4 — 7%, depending upon fiducial region. With the 
addition of central magnetic field and enhanced track- 
ing in the Run II detector, the background is expected 



25 



to be reduced even further. 

To allow direct comparisons of experimental results 
with QCD predictions, the experimental results are 
fully corrected for detector acceptance, experimental 
efficiencies, and detector resolution effects. For the 
Run I data, CDF and D0 used Monte Carlo simu- 
lations to determine the corrections for the detector 
efficiency, acceptance, and resolution as a function of 
Pt and y. CDF generated the di-electron signal using 
PYTHIA [ Hi , processing them with the CDF detector 
simulation and reconstruction programs. Additionally, 
they used PHOTOS [ || to simulate final state QED 
radiation from the J*/Z — * e + e~ vertex. PYTHIA and 
the detector simulation were tuned to obtain satisfac- 
tory agreement with data [ |94). D0 generated 7* jZ 
events using the resummed prediction from LEGACY [ 
[95| , smearing the decay electrons with a parameterized 
detector simulation which included final-state QED ra- 
diation corrections. The parameters and resolutions in 
the detector simulation were tuned to obtain agree- 
ment with data [ |9(| . 

In general, the predictions for vector boson produc- 
tion do not yet include QED effects, therefore the 
experimental corrections must attempt to account for 
them. At a minimum, final state QED radiation must 
be included because the effects are large [ , and CDF 
and D0 included such corrections in their Run I mea- 
surements. As there is yet no numerical implementa- 
tion of the QED corrections analogous to the QCD soft 
gluon rcsummation formalism, initial state QED effects 
have not been considered. Initial state QED radiation 
effects in pt are expected to be similar to those in 
QCD, but since the coupling is much smaller, the pt 
distribution due to QED effects should be much softer. 
It is expected that resummation of the initial state 
soft-photon emission can be implemented similarly to 
the soft gluon case [ |97j . 

Using 110 pb _1 e + e~ data from Run I, both CDF 
and D0 have measured the Drell-Yan cross section 
da/dp T [ 94, [)6) and CDF has measured da/dy [ 
fH ||. For both CDF and D0, the measurement 
error at the peak of the pt distribution (~ 3 GeV/c) 
was ~ 6%. The measurement error at \y\ = was 
~ 5%. The combined efficiency and acceptance for 
these measurements was ~ 33%. 

The dominant systematic uncertainties in these 
measurements are the efficiency and background 
corrections. Generally, the uncertainty in shape is 
more problematic than the uncertainty in overall nor- 
malization. As mentioned, the level of background is 
expected to be reduced in Run II with the enhanced 
tracking available in both detectors. This will also 
reduce the effect of the normalization-uncertainty on 
the final measurement. The uncertainty in the shape of 



the background is dominated by statistics from events 
that satisfy multijet and direct-7 triggers, which also 
satisfy the kinematic and fiducial requirements nec- 
essary for the Z/j* analysis. Again, the enhanced 
statistics of Run II should allow this uncertainty to by 
reduced. The overall electron identification efficiency is 
well-known in both experiments-^ ~ 0.5%. The shape 
as a function of pt is dominated by the isolation re- 
quirement being spoiled by jet activity nearby the elec- 
tron cluster. Understanding this effect requires either 
excellent GEANT-type Monte Carlo or a great number 
of di-electron events from single-electron triggers, so 
one can investigate the effects of hadronic activity on 
electron isolation in an unbiased manner. We expect 
to see improvements in both areas in Run II, hopefully 
reducing the uncertainty in the shape of the efficiency 
as a function of px from ~ 3 — 5% to ~ 1%. 

The expected Run II measurement errors for the 
da/dy and da/dP^ measurement can be estimated 
from the Run I measurement errors by assuming that 
the statistical ("Stat") and systematic ("Syst") errors 
scale as follows: 




where N cv is the number of events in a bin. The scaling 
of the systematic uncertainties should be considered a 
rough parameterization. Scaling the Run I uncertain- 
ties from CDF to an integrated luminosity of 2 fb _1 , 
we obtain predictions for the total measurement un- 
certainty in da/dy (Fig. ^) and da/dpr (Fig. ffHl ). 

With an expected precision of a few percent, the 
Run II data can provide even more stringent con- 
straints on all aspects of our understanding of Z/~f* 
productions, with important implications on our un- 
derstanding of resummation in pQCD, and the preci- 
sion determination of Myj. The precision of the Run II 
da/dy and da/dpT measurements will most likely be 
limited by systematic uncertainties from corrections 
to the detector acceptance/efficiency and resolution. 
These uncertainties will be constrained through the im- 
proved accuracy of Monte Carlo detector simulations, 
and the additional data available in general during 
Run II. 
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3.3. Offstream Searches for Vector Bosons 
3.3.1. Introduction 

Since their discovery in 1983 [ |0(| [l0l[ |l0§ 

and Z bosons have been studied at hadronic colliders 
only via their leptonic decays. As a matter of fact 
the hadronic decays of these particles are generally so 
difficult to separate from the huge QCD backgrounds 
that, after the extraction of a nice mass bump in the 
jet-jet mass distribution by the UA2 collaboration in 
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Figure 44. The expected CDF Run II measurement 
error on da/dy of e + e~ pairs in the mass range 66 - 
116 GeV/c 2 . "Data/Theory" has been arbitrarily set 
to unity. The error is for 2 fb _1 . 
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Figure 45. The expected CDF Run II measurement 
error on da/dp^ of e + e _ pairs in the mass range 66 - 
116 GeV/c 2 . "Data/Theory" has been arbitrarily set 
to unity. The error is for 2 fb _1 . 



1987 [ 104, 105 1, they have laid dormant for quite a 
while as an electroweak physics topic. 

Things have started to change with the increase of 
collider luminosity and dataset size. During Tevatron's 
Run I, hadronic W decays have been successfully used 
by CDF and DO in the discovery and measurement 
of the top quark properties both in the single lepton 
and fully hadronic final states, and a handful of jet-jet 
masses peaking at 80 GcV have been extracted from 
a subset of high-purity tt events (see Fig. [16]). More 
recently, a signal of Z decays to b- quark pairs has 
emerged in the CDF data (Sec. 3.3.3 ). 

With Run II sample sizes it will be possible to 
search for more such hadronic signals, and some of 
them arc expected to start becoming useful tools for 
other physics advances. In fact, their potential as 
calibration tools for the jet energy measurement is 
high, provided that they can be collected by unbiasing 
triggers. This seems particularly likely in the case of 
Z — > bb decays, where the background has been shown 
to be reducible to a manageable size and the hardware 
tools for collecting them with good efficiency and small 
bandwidth concessions are now available; moreover, for 
6-quark jets the absolute energy scale cannot be fixed 
by photon-jet balancing techniques, due to the rarity of 
events with a photon recoiling against a single 6-quark 
jet: a Z peak may then really be our best chance for 
that purpose. Hadronic W decays will also be an ideal 
calibration tool in tt events, but efforts need to be spent 
on finding their signal in independent data samples. 
These may be provided by diboson production pro- 
cesses, where triggering and background issues are less 
problematic. 

3.3.2. Hadronic Decays of W Bosons 

Searches for a W mass signal in inclusive jet triggers 
have been fruitless in Run I data. With respect to 
the lower energy SppS collider, the Tevatron's higher 
center-of-mass energy is a disadvantage for once, be- 
cause in the face of a four-fold increase in signal cross 
section the background from QCD processes increases 
by an order of magnitude, due to its steep behavior 
with respect to parton x. Moreover, no dedicated 
1ow-£t jet triggers were devised either at CDF or D0 
during Run I, given the experiments' focus on the high 
energy frontier; at the very end of the run, however, a 
sector of CDF's central tracking chamber became in- 
operative due to a broken wire, which allowed 1.9 pb _1 
of data to be collected by a high-bandwidth 12 GeV 
dijet trigger. The data thus gathered did not allow 
the extraction of a W peak either, but lends itself to 
fruitful extrapolations to Run II. 

It seems reasonable to investigate the collection of 
events with very low jet Et by special low-luminosity 



27 



I 2.5 



1.5 



0.5 



Shaded: non-tt background 
Dotted: it 175 MC + bg 



50 



100 



150 



200 



250 



300 



Dijet Mass (GeVAr) 



Figure 46. The jet-jet mass distribution of untagged 
jets in the very high purity single lepton tt candidates 
sample collected by CDF requiring two additional b- 
tagged jets. 



runs, which might become an attractive option in the 
event of temporary inoperativeness of tracking detec- 
tors. Here we examine a scenario where CDF II or D0 
gathers some 100 pb _1 of unprescaled data collected 
with a trigger requiring two calorimeter towers above 
4 GeV at Level 1 and two jets with E T > 10 GeV 
at Level 2, back-to-back in phi -about the smallest 
thresholds that do not saturate the bandwidth. About 
two million W/Z — > jj decays could be then collected, 
which would be reduced to 200,000 after optimized 
kinematic cuts requiring two central back-to-back jets 
and little extra-jet activity. Toy Monte Carlo studies 
suggest that in such a scenario the W mass could be 
fit with a ~ 0.5 GeV uncertainty, provided the avail- 
ability of a prescaled sample with looser requirements 
at Level 2: this would yield an understanding of the 
absolute jet energy scale of the detectors to better than 
1%. Such a dataset could then clearly be used also 
for excellent tests of optimization of jet algorithms, 
and thus offer benefits to any search for hadronically 
decaying massive objects. 

3.3.3. Z Decays to 6-Quark Pairs in Run I 

Thanks to the several million Z decays to &-quark 
pairs collected by the LEP I and SLD experiments 
since 1992, the physics of these decays is extremely 



well studied and understood. At a proton-antiproton 
collider that particular process had not been identified 
before, though; therefore the extraction of a signal in 
Run I data was interesting in its own right. Moreover, 
the knowledge of how to extract a Z peak enables a 
careful design of a dedicated trigger for Run II, which 
may allow us to collect a large sample of these events, 
from which the mass distribution can be fit and thus 
insight can be obtained on the absolute energy scale for 
6-quark jets, substantially reducing one of the critical 
sources of systematic uncertainty in the top quark mass 
measurement. 
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Figure 47. Top: the Z — > bb peak in the signal sample; 
bottom: a signal-depleted sample is used to extract the 
background shape. 



A Z — > bb signal was extracted in Run I CDF 
data using kinematic tools and 6-quark vertex tagging. 
Starting from a dataset of about five million events 
enriched of 6-quark decays collected by a single muon 
trigger, a very tight kinematic selection was devised, 
which increased the signal to background ratio by three 
orders of magnitude. 

The main background to the Z decay to b quarks is 
due to direct QCD production of a bb pair, while non-6 
backgrounds yielding muons can be completely elimi- 
nated by requiring the presence of a secondary vertex 
in each of the two jets; this cut also reduces the flavor 
excitation and gluon splitting contributions to b quarks 
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in the sample quite effectively. Most of the direct bb 
pairs are produced at the Tevatron by a gluon fusion 
process, whose high color charge in the initial state 
and color flow topology are distinctive characteristics. 
To exploit the smaller probability of QCD radiation 
in the signal, the two leading jets were required to be 
back-to-back in <fi, and the sum of all other calorime- 
ter clusters in the event was required to be smaller 
than 10 GeV. These cuts selected 588 events, whose 
jet-jet mass distribution was fit to the sum of QCD 
background and Z decay to yield a signal of 91 ± 30 
events, with a S/N ratio at the Z peak equal to 1/3; 
the fit resulted in a Z mass of 90.0 ± 2.4 GeV (see 
Fig. ^) . If the same analysis should be replicated with 
20 times more statistics and no detector improvements, 
this would yield a relative error in the b-jet energy scale 
smaller than 1%. The picture could be even rosier as 
far as statistics goes, due to the extended lepton cover- 
age and improved silicon tracking the Run II detectors 
will be endowed with, but the strong bias due to the 
triggering lepton will make these findings of difficult 
use for generic 6-quark jets. 

Besides its possible use as a calibration tool, how- 
ever, a Z — > bb peak provides a fine testing ground for 
algorithms designed to improve the jet energy measure- 
ment for b jets, which is one of the critical points for 
the discovery of an intermediate mass Higgs boson in 
Run II. The bb final state is the dominant one in Higgs 
boson decay if Mh < 135 GeV/c 2 [EOfl. Our ability to 



extract this particle from the large QCD background 
in Run II (for instance in the Wbb final state, when 
associated WH production is sought) will therefore 
depend critically on the resolution we can attain on 
the Higgs boson mass as reconstructed from the mea- 
sured 6-quark jet energies: both the possibility to see 
a bump in the mass spectrum of jet pairs associated to 
a leptonically decaying W bosons, and the alternative 
option of applying a mass window cut as a selection 
tool for these events, will strictly depend on the actual 
mass resolution. 

The expected resolution for a generic jet-jet reso- 
nance at CDF and DO was roughly a^., = O.lMjj in 
Run I. A relative improvement of this number by 30% 
would significantly extend [ 106 our discovery reach 
for the Higgs boson in Run II. In order to achieve 
that improvement we must study in detail the char- 
acteristics of 6-quark jets emitted in the Higgs decay, 
and use to their utmost the large amount of available 
information provided by the various detector compo- 
nents D0 and CDF II are made up of. For example, 
three-dimensional tracking in the new SVX II detector 
may allow CDF II to infer the momentum of the escap- 
ing neutrinos in semileptonic 6-quark decays, greatly 
improving the energy measurement of the resulting 



jets; this plan will work well in Run II, given the larger 
acceptance for charged leptons from semileptonic de- 
cays provided by the new detectors. Furthermore, the 
possibility of measuring track momenta to higher ra- 
pidity will allow a fruitful use of tracking information 
to improve the calorimetric measurement of jets. 

A detailed study of the observable quantities of b- 
quark jets produced in Z — > bb decays followed by one 
semileptonic b — > /iX decay have been shown to allow 
a sizable reduction of the width of the reconstructed 
bb peak. The quantities found useful for this purpose 
in the CDF analysis were the muon momentum, the 
projection of missing transverse energy along the jet 
axes, and the charged fraction of the jets. 

The muon momentum is needed in the correction of 
jets containing a semileptonic decay of b quarks, be- 
cause the minimum ionizing muons do not contribute 
linearly to the energy measured in the calorimeter. 
The missing E T , projected along the jet directions in 
the transverse plane, provides useful information on 
the amount of momentum taken away by the neutrino 
in the muon jet and on possible fluctuations of the 
energy measurement in both the muon and the away 
jet. The charged fraction of the jets, defined as the 
ratio between the total momentum of charged tracks 
belonging to a jet and the energy measured in the 
calorimeter, also helps reducing the uncertainty in the 
energy measurement. 

By properly accounting for the value of these observ- 
ables, it was possible to reduce the relative uncertainty 
in the dijet mass measurement, afci/Mjj, by nearly 
50% (see Fig. ||). 

If the alternative plan (described in the following 
section) based on collecting Z — ► bb decays by trig- 
gering directly on secondary vertices in jets at Level 2 
should fail, CDF II will anyways be able to observe a 
peak of several thousand events in the inclusive lep- 
ton datasets by simply replicating the Run I analysis. 
These peaks will not allow a precise calibration tool for 
inclusive 6-quark jets, given the biasing semileptonic 
decay of one of the two jets, but they will anyways be 
extremely useful as a testing ground for the algorithms 
now under development which aim at reducing the jet 
energy resolution. 

3.3.4. Z Decays to &-Quark Pairs in Run II 

A strategy to collect Z — ► bb events in Run II has 
been studied having in mind the CDF trigger config- 
uration. We tried to keep the trigger requirements 
as simple as possible and optimize them in order to 
have acceptable trigger rates (less than 2% of the total 
bandwidth at each level) and maintain a good signal 
efficiency. 

Most Z — > bb events should contain two recon- 
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PYTHIA Z -> bb: Mass Reconstruction 
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The computing power of the Level 3 processors 
should allow the reconstruction of secondary vertices 
online. We require two jets, each with uncorrected 
Et > 10 GeV in a cone of radius R = 0.7, and check 
that the jets contain two displaced vertices. 

A summary of signal efficiency, trigger cross section 
and trigger rate at the three trigger levels is reported 
in Table |] for a typical instantaneous luminosity C = 
1.4- 10 32 cm- 2 s- 1 . 



Table 4 

Summary of the Z — ► bb trigger efficiency, cross section, 
and rate at each level. 



Trigger 


es (%) 




R T (Hz) 


Level 1 


17.0 ±0.2 


(5.7 ±0.3) fib 


800 ± 42 


Level 2 


3.10 ±0.03 


(32 ± 3) nb 


4.4 ±0.4 


Level 3 


2.4±0.1 


(3.2 ±0.9) nb 


0.5 ± 0.1 



Figure 48. The four Gaussian fits show the improve- 
ment of the mass reconstruction for simulated Z — > bb 
events (Pythia V5.7, CDF detector simulation) when 
the observable characteristics of the 6-quark decays are 
properly taken into account in the mass reconstruction. 



Since a(pp -> Z -> bbX) ~ 1 nb, about 2 • 10 6 
Z — > bb events will be produced in Run II, assuming 
2 fb _1 of data collected in two years. We expect to 
collect 48,000 signal events in a total of about 6.4 • 10 6 
with this trigger, corresponding to a S/B of 0.0075 and 
a significance S/VB ~ 19. 



structable secondary vertices in the final state. In 
Run II it will be possible to trigger on secondary ver- 
tex information. Exploiting this feature, it should be 
possible to collect a sample of Z — > bb events without 
particular requirements on the b decay. One purpose 
of this dataset is to provide a sample to calibrate the 
calorimeter energy scale for jets containing b quarks, 
which may limit our ability to measure the top mass. 
Therefore the trigger path should bias the energy scale 
measurement as little as possible. 

At Level 1 the XFT (eXtremely Fast Tracker) in- 
formation will be available: transverse momentum Pt, 
azimuthal angle 4> and charge sign of particles crossing 
all of the Central Outer Tracker (COT) layers. We 
require the presence of two central high Pt tracks in 
opposite hemispheres. Cuts are chosen on the track 
Pt and A<f) to maximize the statistical significance 
(S/y/B). The track P T cuts are 6 GeV and 4 GeV, 
with A4> > 150°. 

At Level 2 the rate is reduced by requiring the tracks 
to have a finite impact parameter. The SVT (Silicon 
Vertex Tracker) processor will provide this informa- 
tion. The best significance is found by requiring two 
SVT tracks with 120 \xm < \d\ < 1000 /im. 
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Figure 49. Mjj distribution expected from Run II: 
signal±background. 



It is also advantageous to require the presence of no 
additional jets with Et > 10 GeV in the event, either 
at trigger level or at a pre-analysis level. This would 
increase the signal fraction to 0.013 and the significance 
to S/ \[B = 24. We simulated the signal extraction 
procedure under these assumptions. We searched for 
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the Z — > bb signal in the jet-jet invariant mass distri- 
bution Mjj. We expect to see an enhancement cor- 
responding to the Z mass, since the statistical signifi- 
cance of the signal is high. Outside of the signal region 
the Mjj spectrum should be approximately described 
by a decreasing exponential. 

The PYTHIA Monte Carlo is used to model the 
Z — ► bb invariant mass spectrum. The background Mjj 
distribution expected for the background is inferred 
from generic Run I dijet data by assuming that the 
request of secondary vertices does not modify its shape. 
The Mjj distributions for signal and background nor- 
malized to 2 fb^ 1 have been added (see Fig. |49| ). 
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Figure 50. Background-subtracted Mjj distribution 
for Z events obtained from a pseudo-experiment. 



Because of the very low S/B ratio, we cannot dis- 
tinguish by eye the signal peak in Fig. We con- 
ducted 1000 pseudo-experiments, varying every time 
the bin content according to a poissonian distribu- 
tion with mean equal to the bin content itself. For 
each pseudo-experiment the region outside the Z peak 
was fitted with a decreasing exponential. The value 
of this function was subtracted, bin by bin, from the 
total spectrum. The results obtained in one of these 
pseudo-experiment subtractions is shown in Fig. |5^, 
together with a gaussian fit to the excess found. With 
this method we typically obtain for the jet-jet invari- 
ant mass peak position M meas = 84.0 GeV/c 2 , with a 
statistical uncertainty of 1.5 GeV/c 2 on Af meas .^| 

The interval of ±2er around the Z peak 
(55 GeV/c 2 < Mjj < 115 GeV/c 2 ) contains about 
93% of the signal events and 25% of the background 

^ The jet energy corrections applied here are not optimized for 
b jets, whose higher mass and decay properties are sensibly 
different from those of generic light-quark and gluon jets: the 
invariant mass peak is thus about 7 GeV lower than the true Z 



events. Therefore in this region S/B = 0.05 and the 
statistical significance is S/V~B — 44. 

The error on the Z mass can be directly translated 
into an error on the 6-jet energy scale. Measuring the 
Z mass with an error of the order of 1.5 GeV/c 2 will 
allow to determine the 6— jet absolute energy scale with 
an uncertainty of 1.7%. 

3.3.5. Other Hadronic Signals 

Besides its interest for the study of trilinear gauge 
boson couplings, discussed in Sec. |[ associated pro- 
duction of two vector bosons yielding a leptonic and 
a hadronic decay may provide additional handles for 
the physics of hadronic resonances, both because of 
the ease of collecting these events with good efficiency 
in high- Ft lepton triggers and because of the larger 
signal to background ratio with respect to single boson 
production processes, due mainly to the reduction of 
background processes with gluons in the initial state. 

WW production is the best example: in a sample 
of 5 fb -1 , for instance, about 3000 pp — > WW — ► Ivjj 
events can be collected by applying standard cuts on 
the leptonic decay products and requiring two jets with 
uncorrected Et > 15 GeV; the signal to noise ratio 
is then close to 1/40 before any optimized selection. 
Thence an observation of the W peak in the dijet mass 
distribution will be relatively easy to obtain. System- 
atic effects in the mass fits due to the low S/N ratio 
may make this sample of little impact as a source of 
knowledge of jet energy scale when compared to the 
high-purity W — > jj samples that single lepton tt de- 
cays may provide; but the signal may still be extremely 
useful for the study of jet resolutions. 

Another process that will be likely observed in 
Run II is associated WZ production with a subsequent 
leptonic decay of the W boson and a decay of the 
Z to 6-quark jets. In 5 fb -1 about 500 such events 
can be collected by the lepton triggers, from where 
secondary vertex tagging can considerably increase the 
signal purity. Despite its small size, this signal may be 
of fundamental importance in checking systematic un- 
certainties in the standard model Higgs boson search, 
which will mainly focus on the very same dataset with 
very similar analysis cuts. 

Finally, the possibility of collecting Z — > bb in pho- 
ton triggers has been investigated. Zj production has 
a small cross section -about 2.5 pb for a photon with 
Et > 10 GeV and two &-quark jets (see Fig. jsiT)- 
but the heavy flavor decay provides a quite distinc- 
tive signature, and the process could be easily put in 
evidence if a sufficiently low trigger threshold were set 
on the photon transverse energy, or alternatively if a 76 
trigger could be devised. The advantage of the process 
producing an additional photon over inclusive Z pro- 
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Cross Section vs Pt 
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Figure 51. The cross section times branching ratio 
for Z7 production with subsequent Z —* bb decay, for 
different choices of PDF sets, c omp uted with the Baur- 
Berger-Zeppenfeld generator 
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By probing the spin structure of W production, this 
measurement provides another method that is indepen- 
dent of purely QCD analyses, while adding to the list of 
measurements using vector bosons to study NLO QCD. 
Moreover, the measurement of the angular distribution 
of the decay leptons is also of importance for the W 
mass measurements, because the next-to-leading or- 
der QCD corrections to the angular distribution are a 
non-negligible contribution to the W mass. 

The measurement of the angular distribution of elec- 
trons from W b osons ob tained with Run I data col- 
lected by D0 [ |l!0[ 111 ] is statistically limited (see 
Fig. ^||). While a calculation that includes QCD effects 
is preferred over one that does not, this preference is 
not strong enough to exclude a pt independent angu- 
lar parameter oli- With the next collider run starting 
in the near future, it is worthwhile looking at the sen- 
sitivity of this measurement in Run II. In the following 
discussion we will estimate the size of statistical and 
systematic errors to this measurement with 2 fb _1 in 
Run II. 



duction lays in the suppression of the irreducible back- 
ground from direct QCD bb production: since initial 
state photons can only be present in quark-antiquark 
annihilation diagrams, at Tevatron energies the bbj 
background is suppressed by roughly an order of mag- 
nitude more than the Z signal with respect to the 
searches in inclusive bb samples. It is however difficult 
to devise an unprescaled trigger capable of a sufficient 
efficiency for the signal while maintaining the total rate 
at an acceptable level. In particular, a 76 trigger would 
have to require at Level 1 both a photon candidate and 
one or two charged tracks; the tracks and the photon 
would then need to be separated in azimuthal angle, 
to reduce the rate of fake photon signals from QCD 
events. There is currently no plan to devise such a 
trigger in CDF II or D0 , while the lowest unprescaled 
photon triggers will collect events with Et 1 > 25 GeV, 
where the 766 cross section is only 0.6 pb. With 5 fb _1 
it will be relatively easy to isolate a signal of one or two 
hundred events over a similarly sized background, but, 
given the small size, its exploitation appears dubious. 
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3.4. Lepton Angular Distributions in W Boson 
Decay 

Next-to-leading order perturbative QCD predicts 
that in W — > Iv decays a n an gular distribution of 
(1 ± aicos0* + a 2 cos 2 #*) [ |l08| should be observed, 
where 9* is the polar angle of the decay lepton in the 



Collins-Soper frame [ 109 . In the presence of QCD 



corrections, the parameters ol\ and a 2 are functions of 
p^ , the W boson transverse momentum. 

The measurement of oti serves as a probe of NLO 
QCD, using the well understood I^-fermion coupling. 



Figure 52. D0 Run I result: the measured as a 
function of px with combined statistical and system- 
atic errors compared to the NLO QCD calculation by 
Mirkes (curve) and calculation in the absence of QCD 
effects (horizontal line). The statistical errors alone 
are shown as horizontal ticks. 



The expected statistical errors in the determination 
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of «2 should simply scale like the inverse of the square 
root of the number of events. We consequently need 
to calculate the expected number of W boson events 
under Run II conditions. We get a factor of 57 in W 
boson statistics which breaks down as follows (see [ 
0): 



Nw,RunII , j, N e + 

- =Jlum X J"2.0TeV X 



W,RunI 



N,, 



^tracking 
€no tracking 



(3) 



where fi urn — 20 indicates the increase in luminosity, 
/2.0TeV = 1.2 indicates the increase in W cross section 
due to the increase of the center-of-mass energy from 
y/s = 1.8 TeV to 2 TeV, (N e + N^)/N e = 2 is the addi- 
tional statistics gained by including the muon channel, 
and 



^tracking 
€no tracking 



0.95 
~08~ 



(4) 



is the increase in efficiency due to tracking capabilities 
of the upgraded D0 detector. 

The statistical errors of the Run I measurement are 
therefore scaled by 1/V57, as shown in Fig. The 



is therefore necessary. Table [| shows a summary of 
statistical and systematic errors for Run I and Run II. 
Since the modeling of the hadronic recoil is done from 
Z data, the error due to the hadronic resolution will 
scale with Z statistics up to a point. The estimate of 
this error is done by scaling the number of Z events 
by a factor of y/28.5 (the same as for W events but 
excluding muons). The error on the hadronic response 
will also improve with increased Z statistics. Conse- 
quently, the largest errors left are the ones due to the 
QCD multi-jet background and the electromagnetic 
scale. It will be of crucial importance for this and other 
electroweak measurements to better estimate the QCD 
background. This goal can partially be reache d b y 
taking more events with a QCD monitor trigger [ 111 ] . 
Of course, it would be best if the QCD background 
fraction could be reduced even further. At this point 
it is not obvious if this is feasible since the current set 
of electron identification cuts is already very efficient 
in reducing this background. A Monte Carlo program 
that correctly models QCD multi-jet events including 
a realistic detector model could also help in better es- 
timating the shape of this background. For the low pt 
region, the dominant error will be the electromagnetic 
scale. There is however currently no good estimate by 
how much this error will be reduced. 



no QCD effects 
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Figure 53. Estimated sensitivity of the a 2 measure- 
ment obtained by scaling statistical errors to Run II 
conditions. 



statistical uncertainties to this measurement become 
quite small, and a further look at systematic errors 



Table 5 

Central values and statistical errors for a.2 and sys- 
tematic errors due to backgrounds, the hadronic en- 
ergy scale and resolution, and the choice of ct\. All 
systematic errors are for Run I if not otherwise noted. 



p T [GeV] 


0-10 


10 - 20 


20 - 35 


35 - 200 


mean px [GeV] 


5.3 


13.3 


25.7 


52.9 


c*2, Run I 


1.07 


0.82 


0.49 


0.10 




±0.13 


±0.25 


±0.37 


±0.37 


c*2, Run II 


0.98 


0.89 


0.68 


0.24 


Stat, errors, e 


±0.024 


±0.047 


±0.069 


±0.069 




±0.017 


±0.033 


±0.049 


±0.049 


total syst error 


±0.08 


±0.09 


±0.12 


±0.12 


QCD 


±0.04 


±0.05 


±0.09 


±0.07 


Z 


±0.01 


±0.02 


±0.02 


±0.04 


top 


±0.00 


±0.00 


±0.00 


±0.02 


EM scale 


±0.06 


±0.05 


±0.03 


±0.04 


had scale 


±0.03 


±0.01 


±0.04 


±0.04 


had. resolution 


±0.02 


±0.02 


±0.05 


±0.06 


had. res. Run II 


±0.004 


±0.004 


±0.009 


±0.011 


<X\ 


±0.01 


±0.05 


±0.03 


±0.03 



In the current measurement we had to fix a± to the 
value predicted by the QCD calculation since even after 
summing over both W signs we are slightly sensitive 
to ai due to acceptance effects. The error due to this 
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choice for ot\ was estimated by setting a,\ = 2.0, the 
value predicted by the V — A theory in the absence of 
QCD effects (see Table §. This error is non-negligible. 
Since the central magnet in Run II will allow for sign 
identification of electrons, u\ and U2 could thus be 
measured simultaneously, eliminating the need for the 
above assumption for a%. This will reduce the error 
due to the choice of ot\. While this is a nice extension 
of this measurement, it is not clear at this point by 
how much it will numerically improve the significance 
of the measurement of a.2- 

In the above estimate of the errors, the binning in 
Pjf used for the Run I measurement was kept un- 
changed. With larger statistics, one would probably 
choose a finer binning in p^ , allowing for bins with 
larger mean p^ . This would increase the sensitivity in 
the area where the deviation of the angular distribution 
due to QCD effects is most pronounced. 

The Monte Carlo used in the current analysis, which 
was originally developed for the measurement of the W 
mass [ |8) at D0, treats hadronic jets as point particles 
and the hadronic recoil is treated as a single jet. This 
is clearly a simplification of the true processes involved 
and a real next-to-leading order event generator would 
be useful. 

In addition to the experimental improvements dis- 
cussed thus far, this measurement will be sensitive to 
W production models. These models have to be con- 
strained by independent measurements. 

To summarize, in Run II the measurement of the 
angular distribution of electrons from W boson de- 
cays will be systematically limited. While the recoil 
response and resolution will improve with increased Z 
statistics, the estimate of the QCD background frac- 
tion and shape becomes a limiting factor. It is difficult 
to estimate by how much the other dominant error, 
the error due to the uncertainty of the electromagnetic 
scale, will be reduced in Run II. Other improvements 
not quantified here are expected from a finer binning 
in pt and sign identification of electrons. 



plings. Deviations of the couplings from the SM values 
would indicate the presence of new physical phenom- 
ena. 

The purpose of this section is to present a brief 
overview of recent theoretical advances in understand- 
ing di-boson production in hadronic collisions, and to 
highlight Run II opportunities for studying the physics 
of vector boson pair production. Because of the large 
anticipated size of the data sample, JCdt = 2 fb _1 , 
interesting processes and final states that were not 
studied in Run I will become available. These are 
discussed, as well as prospects available in extensions 
of the Run I Wj, WW, and WZ analyses to Run II. 
This is meant to be an improvement over the forecasts 
of the TeV_2000 Report [ [r|, which was written in 
1995 before we had the benefit of having performed 
the Run lb analyses. Indeed, some of the TeV_2000 Re- 
port's prognostications for Run II limits were achieved 
in Run I. 

We begin with a brief summary of the trilinear 
gauge boson couplings and how they are parameter- 
ized. Next, we give a short description of new theo- 
retical developments. Following that, we summarize 
the anomalous coupling limits obtained in Run la and 
Run lb, and compare the Run lb results with what 
we expected we would obtain, based on a simple ex- 
trapolation from Run la. This exercise in hindsight 
provides both a calibration for, and a cross-check of, 
the extrapolation method. The subsequent section 
provides expectations for anomalous coupling limits 
from the Run II analyses based on extrapolation of the 
Run I analyses to higher integrated luminosity. Next, 
we provide comments on some Run II analyses, and, 
lastly, discuss new channels and analyses which will 
become feasible in Run II, in particular the prospects 
for measuring the ZZV couplings via ZZ production. 



4. Vector Boson Pair Production and Trilin- 
ear Gauge Boson Couplings Prospects for 
Run II§ 

The Standard Model of electroweak interactions 
makes precise predictions for the couplings between 
gauge bosons due to the non-abelian gauge symme- 
try of SU(2)l (g) U(1)y- These self-interactions are 
described by the triple gauge boson (trilinear) WWj, 
WWZ, Zjj, ZZj and ZZZ couplings and the quartic 
couplings. Vector boson pair production provides a 
sensitive ground for direct tests of the trilinear cou- 



4.1. Trilinear Couplings 

The WWV (V — 7 or Z) vertices ar e described by a 
general effective Lagrangian [ 112j , 
couplings, gww-y = — e and gwwz 
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with two overall 
and 



§ Contributed by: U. Baur, H.T. Diehl and D. Rainwater 



-e ■ cot 9w 

six dimensionless couplings g\ , Kv, and Ay (V = 7 
or Z), after imposing C, P, and CP invariance. The 
Wa (q) Wt (q) Vfj, (p) vertex function (where all mo- 
menta are outgoing, p + q + q = 0) in presence of 
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non-standard couplings is given by: 



factors for these couplings are 



^«/3f / 
WWV -a 8u I V i i \ 

gwwv \ 



M 2 W 



«V" (gY + Kv + Xv 



Here, M^y is the VF-boson mass. Electromagnetic 
gauge invariance requires that gj — 1, which we as- 
sume throughout this paper. The SM Lagrangian is 
obtained by setting gj = gf = 1, Ky = 1 (A/ty = 
Ky — 1 = 0) and Ay = 0. 

A different set of parameters, motivated by SU(2) x 
U(l) gauge inv arian ce, had been used by the LEP 
collaborations [ |l!4| prior to 1998. This set consists 
of three independent couplings otBcj,, ccw4> an d aw'- 
OiB4> = Ak 7 — Agf cos 2 9\y, aw<j> = Agf cos 2 &w and 
aw = A-y. The remaining WWZ coupling parameters 
Xz and Akz are determined by the relations Xz = A 7 
and Ak z = -Ak 7 tan 2 6 W + A</f. The HISZ rela- 
tions [115] which have been used by the D0 and CDF 



collaborations are also based on this set with the ad- 
ditional constraint olb^> — awtj>- 

The di-boson production cross sections with non-SM 
couplings grow with the parton center of mass energy 
V§' In order to avoid violation of S'-matrix unitarity, 
the anomalous couplings a = gY , A«y, Ay are taken 
as momentum dependent form factors with a scale App 



a(s) 



(1 



s/^ppY 



(G) 



and n — 2 (dipole form factor). 

The Z a (q 1 ) 7 /3 (g 2 ) V»(P) (V = 7 or Z) vertices con- 
tributing to Z7 production are described by a general 
vertex function [ 112| with eight dimcnsionlcss cou- 



plings hY (i = 1, . . . , 4 ; V = 7 or Z): 



Z-iV — 



Ml 

h%_ 

+ M% 



h\ ( q $g af3 - q %g^) 

P a ((P-q 2 )g^-q2pP) 



(7) 



Q2p 



In the SM, all hY's are zero. The couplings h\ and 
hX violate CP; all couplings are C-odd. The form 



hi 



{l + s/A% F Y 



where one usually assumes that n 
n = 4 for i 



(8) 



3 for i = l,3 and 



2,4 [ |107|. 

In the SM, the £ + £~j final state can be produced via 
radiative decays of the Z boson or by production of a 
boson pair via t- or it-channel quark exchange. The 
former process is the dominant source of events with 
small opening angle between the photon and charged 
lepton and for events with a low value of photon 
transverse energy, Ej,. Events produced by the latter 
process have lepton-pair invariant mass, mu, close to 
Mz and three-body invariant mass, m« 7 , larger than 
Mz- Anomalous ZZj or Z77 couplings would enhance 
the cross section for Zj production, particularly for 
high-i?y photons, relative to the SM expectations. 

The most general form of the Z ct (q 1 ) Z^(g 2 ) V»{P) 
vertex function can be written in the form [ pL 12| 



ZZV 



P 1 



Ml 



+ if< 



V uaf3p 



(qi - 92), 



(9) 



CP invariance forbids fY and parity conservation re- 



quires that f¥ vanishes. 



In the SM, f\ 



0. 



.S'-matrix unitarity requires a form factor behavior for 
ZZV couplings similar to that of h\ and h\ [ 116 



(Eq. (|) with n = 3) 

Although the WWV, ZjV and ZZV couplings usu- 
ally are assumed to be real, they are in general complex 
quantities. 

In theories which go beyond the SM, the WWV cou- 
plings are expected to be at most 0(M W / A 2 ) where 
A ~ App is the scale of new physics. Z7V and ZZV 
couplings are at most 0(M Z / A 4 ). 

4.2. Recent Theoretical Developments 
4.2.1. Pa rame terization of ZjV Couplings 

In Ref. [ 117 it was pointed out that the couplings 
hY have to be purely imaginary quantities in order to 
guarantee that an effective Lagrangian which would 
lead to a vertex function of the form of Eq. (^) is 
hermitian. In contrast, the WWV and ZZV couplings 
are normalized such that real couplings automatically 
correspond to a hermitian effective Lagrangian. Since 
one usually assumes real couplings when placing limits 
on anomalous vector boson self-couplings, it is useful 
to replace the couplings hY by 



for the Run II analyses and beyond. 



(10) 
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Figure 54. The photon transverse momentum distri- 
bution in pp — > £ + £~^ at the Tevatron in the SM and 
for anomalous ZZj couplings. The cuts imposed are 
described in the text. 



h't 's and /i^'s of equal magnitude result in virtually 
the same differential cross sections at high energies. 
This is illustrated in Fig. |34| for the case |/i'f | = 
\hf \ — 0.3. In order to simulate detector response, 
the following cuts have been imposed in Fig. 



Pr(j) > 10 GeV, 
p T (£) > 20 GeV, 
m(U) > 75 GeV, 

and 

AR(£-y) > 0.7. 



|J7(7)I<2.5, (11) 
< 2-5, (12) 
m(lli) > 100 GeV, (13) 



(14) 



The form factor scale has been chosen to be Aff = 
750 GeV. 

Unlike for real h^ 4 couplings, the interference terms 
between the SM and the non-SM contributions do not 
vanish in the squared matrix element for real values 
of h' 3 4 . Thus, for intermediate values of Pt(i), the 
differential cross sections for values of h' 34 of equal 
magnitude but opposite sign slightly differ. Since most 
of the sensitivity to anomalous couplings originates 
from the high energy domain, the limits for h'^f are 
expected to be almost identical to those obtained for 
hY . In the following we therefore list limits only for 

4.2.2. NLO QCD Corrections to Vector Boson 
Pair Production 

In the Run I di-boson analyses, data were compared 
with leading order production calculations to extract 
limits on the WWV and Z^V couplings. The effect 



Figure 55. The e + e~ transverse momentum distri- 
bution in pp — > W + W~ + X — > e + e~p' T + X at 
the Tevatron in the SM. The following cuts were im- 
posed to simulate detector response: pr(e) > 20 GeV, 
\f](e)\ < 2.5 and i> T > 20 GeV. For the NLO 0-jet 
curve (dotted line), jets with pr(j) > 20 GeV and 
\r](j)\ < 3.5 were vetoed. 



of higher order QCD corrections was simulated by 
multiplying the differential cross sections by a simple 
constant fc-factor 



k = l + —a s 



(15) 



NLO calculations have shown [ |ll§[ that the 0(a s ) 
QCD corrections in the SM depend logarithmically on 
s and become large at high energies, due to gluon- 
induced partonic subprocesses, which only enter at 
NLO. An example is shown in Fig. |5^, where we display 
the transverse momentum distribution of the charged 
lepton pair \ripp -> W + W~ +X -> e+e~j4 T + V. NLO 
corrections are seen to be very large (0(10)) at high pr, 
and dramatically alter the shape of the distribution. 
Qualitatively this is precisely what one expects from 
non-standard WWV couplings. Since the real emission 
diagrams are responsible for the increase of the QCD 
corrections with pt, a jet veto drastically reduces the 
size of the QCD corrections (dotted line). It should 
be noted that the NLO QCD corrections reduce the 
W + W~ cross section when a jet veto is imposed. 

With the more than 20-fold increase in statistics ex- 
pected in Run II, it is clear that the QCD corrections 
to di-boson production must properly be taken into 
account when information on anomalous couplings is 
extracted. Over the pas t seven yea rs the N LO QCD 
corrections to Wj \ |U9|, Z7 [[l20), WW [ [ml and 



WZ production [ 122 including non-standard WWV 



3G 



Table 6 
The pp - 



w+w- 



X 



X cross section 



for y/s = 2 TeV with p T (e) > 20 GeV, |?7(e)| < 2.5, 
and p T > 20 GeV. Jets are required to have pr(j) > 
20 GeV, \rj(j)\ < 2.5. Results are shown for the calcu- 
lations of Ref. [ |m| (BHO) and Ref. [ ft24j (DKS) 
with factorization scale Q 2 = M^, and using the 
CTEQ4M [ 126] set of parton distribution functions. 





BHO 


DKS 




a (fb) 


a (fb) 


Standard Model 


Born 


61.2 


61.2 


NLO 


80.9 


81.3 


NLO 0-jet 


65.6 


65.3 


A ff f = 0.5, X z 


= A 7 = 0.1, An z 


= Ak 7 = 0.3 




A FF = 2 TeV 




Born 


82.7 


82.8 


NLO 


106.5 


107.0 


NLO 0-jet 


84.2 


83.7 



and ZjV couplings have been calculated using the 
narrow width approximation and ignoring spin corre- 
lations in the finite virtual corrections. Recently, more 
complete calculations have become available which 
properly take into a ccou nt the previously ignored spin 
correlations. 



Ref. 



123] also includes single resonant 
diagrams and finite W and Z width effects in the cal- 
culation; however, no anomalous couplings are taken 
into account. Refs. [124, 125 use the narrow width ap- 
proximation, but do include the option of non-standard 
WWV couplings. 

The contribution of the finite virtual corrections to 
the NLO cross section is smaller than about 10% for 
all di-boson processes. The spin correlations ignored 
in Refs. [ |TT|, |l2^, [l2l], ^ therefore are expected to 
have a rather small effect on the total cross section, 
as well as on most distributions. This expectation is 
confirmed by an ex plici t comparison betw een the cal- 
culations of Refs. [|l24| (DKS) and [ |p (BHO) for 
W + W~ production which is shown in Table |[ The 
two calculations are seen to agree at the 1% level. 

A more detailed and careful comparison between the 
BHO and DKS calculations for W+W~ and W ± Z 
prod uction at the LHC has been carried out in Ref. [ 
|127| . The WZ cross sections of the two calculations 
were found to agree within 1.5%, whereas in the WW 
case deviations of up to 3.8% were observed at NLO. 
Subsequently, a s mall error in the BHO W + W~ code 
was discovered 



128 



After correction of this error, 
the BHO and DKS calculations of W + W~ production 
at the LHC agree to better than 0.5% [ p8|. 



4.3. Summary of Run la and Run lb Analyses 
at CDF and D0 

This subsection contains a summary of the published 
CDF and D0 analyses. The intent is to provide an 
overview of the progress, a list of published papers, 
and tables which provide the basis for the comparison 
made in the next subsection. 

4.3.1. WWj and WWZ Couplings 

The D0 and CDF collaborations have performed 
several searches for anomalous WWj and WWZ cou- 
plings. Studies [ pM EM p| of pp -> Wj + X 



have shown that the transverse energy spectrum of 
the photons agreed with that expected from SM pro- 
duction. Searches [ 132, 133, 1 34 1 for an excess of 
pp — » WW + X, where the W bosons each decayed to 
Iv {I = e or fi) , yielded events which matched the SM 
prediction. Furt her, the px spectrum of the charged 
leptons agreed ] |l 34f| with the prediction. Studies [ 
1351 , HI [137|, |138| of the processes pp -> WW + X and 
pp — > WZ + X, where one W boson decayed to a lep- 
ton or anti-lepton and the corresponding anti-neutrino 
or neutrino and the other vector boson decayed to a 
quark-antiquark pair manifested as jets, yielded no ex- 
cess of events and a W boson transverse energy spec- 
trum which matched the expected background plus 
SM signal. Lastly, D0 studied [ 135] the process 
pp — > WZ + X where the Z boson decayed to ee and 
the W boson decayed to either ev or [iv. Limits on 
anomalous WW^f and WWZ couplings were derived 
from each of these analyses. Several [ 130, 132, 136 
of these analyses were presented in detail in Ref. [ 
13£]. The results of all of the D0 analyses were com- 
bined [ [L 381, using the maximum-likelihood method [ 



1391, p-401, to form the Tevatron's most restrictive limits 
on anomalous WW~f and WWZ couplings. 

Table [7] shows the anomalous coupling limits 
achieved in each of the analyses described above, the 
luminosity used, and the reference to the paper in 
which the result was published. It should be noted 
that many of the papers published limits under sev- 
eral assumptions for the relations between the coupling 
parameters and with several values of the form factor 
Aff- Only those limits from the case A = A 7 = Xz and 
Ak = An 1 = Akz are listed, except for Wj and WZ 
where only WW^f and WWZ couplings, respectively, 
are relevant. 

4.3.2. ZjV Couplings 

The D0 and CDF collaborations have also per- 
formed several searches for anomalous Zjj and ZZ~/ 
couplings. 



Studies 

x -> e+t 



141 



142, 143 | of the process 



pp — > Zj + X — > £ + £~j + X have shown that the event 
yield and transverse energy spectrum of the photons 
agreed with that expected from SM Zj production, 
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Table 7 

95% confidence level WW^f and WWZ anomalous coupling limits achieved in Run I analyses by the D0 and CDF 
Collaborations. 



Analysis 



Ref. 



Run 



Lum. (pb 1 ) 



A 



FF 



A.C. Limit (95% CL) 



CDF Wj -> 
evj and fii/j 



la 



20 



1.5 TeV 



-0.7 < A < 0.7 
-2.2 < Ak < 2.3 



-0.6 < A < 0.6 
-1.6 < Ak < 1.8 



D0 W-/^ 
ev"f and fivj 



130 1 



la 



13- 



1.5 TeV 



D0 W 1 -> 
evj and [ii/j 



la + lb 



92.8 



1.5 TeV 



-0.31 < A < 0.29 
-0.93 < Ak < 0.94 



D0 WW 
Dilepton 



132 1 



la 



14 



900 GeV 
Equal Couplings 



-2.1 < A < 2.1 
-2.6 < Ak < 2.8 



CDF WW 
Dilepton 



133 



la + lb 



108 



1.0 TeV 
Equal Couplings 



-0.9 < A < 0.9 
-1.0 < Ak < 1.3 



D0 WW 
Dilepton 



la + lb 



97 



1.5 TeV 
Equal Coupling 



-0.53 < A < 0.56 
-0.62 < Ak < 0.77 



CDF WW and WZ 
leptons + jets 



135 1 



la 



19.6 



1.0 TeV 
Equal Couplings 



-0.81 < A < 0.84 
-1.11 < Ak < 1.27 



D0 WW and WZ 



136 



la 



13.7 



1.5 TeV 
Equal Couplings 



-0.6 < A < 0.7 
-0.9 < Ak < 1.1 



D0 WW and WZ 



137| 



la + lb 



96 



1.5 TeV 
Equal Couplings 



-0.36 < A < 0.39 
-0.47 < Ak < 0.63 



D0 WW and WZ 



137 



la + lb 



96 



2.0 TeV 
Equal Couplings 



-0.33 < A < 0.36 
-0.43 < Ak < 0.59 



D0 WW and WZ 



138 1 



lb 



81 



2.0 TeV 
Equal Couplings 



-0.43 < A < 0.44 
-0.60 < Ak < 0.74 



D0 WZ — > 
ee\iv and eeev 



138 



lb 



92 



1.0 TeV 



-1.42 < X z < 1.42 
-1.63 < A ff f < 1.63 



D0 Combined 



138 1 



la + lb 



96 



2.0 TeV 
Equal Couplings 



-0.18 < A < 0.19 
-0.25 < Ak < 0.39 



though it is noted that there were two Zj — > ee-j events 
with photons of Et greater than 70 GeV, expected in 
only 7.3% of trial experiments. 

Studies of the process Zj — > vDj have the ad- 
vantage of the higher branching fraction for decay to 
neutrinos than does the charged-lepton decay mode. 
Furthermore, there is no final state radiation because 
the neutrinos are electrically neutral. However, the 
signal-to-background ratio is rather lower than in the 
charged-lepton analysis. D0 has published [ 144 the 
results of the Run la Zj — > vv^( analysis. Again, the 
spectrum of the transverse energy of the photons, for 
E% > 40 GeV, agreed with the SM prediction. 

D0 and CDF produced limits on anomalous ZjV 
(V = Z."f) couplings using a fit to the spectrum. 
The D0 Run la and Run lb results were combined 
in Ref. [ 143 1. Table shows a compilation of all the 



Run I CDF and D0 results. 



The limits for /ig (hj ) 



and hYo (h^o) are almost identical to those obtained 



for /ig (/ifo) an d are > therefore, not shown. 



4.4. Hindsight: Extrapolating Run la Results 
to Run lb 

It is interesting to see how well one can "predict" 
the Run lb limits based on the Run la results and a 
simple rule for scaling the limits based on the increase 
in the luminosity. 

The WW~i and WWZ anomalous coupling limits 
should scale by 



Cdt 



1/4 



(16) 



One square-root comes from the decrease in the statis- 
tical uncertainty of the cross section (as a function of 
Et, for instance) and the other from the fact that the 
differential cross section is a quadratic function of the 
anomalous couplings. 

The Z77 and ZZ^f anomalous coupling limits would 
also scale by the fourth-root of the ratio of the lumi- 
nosities, except that the limits depend very strongly 
on the form-factor scale. 
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Table 8 



95% confidence level ZZ^ and Z77 anomalous coupling limits achieved in the Run I analyses by the D0 and CDF 
Collaborations. 



Analysis Ref. 


Run 


Lum. (pb 1 ) 




A.C. Limit (95% CL) 


CDF Z7 -> [ 


141 


la 


20 


500 GeV 


-3.0 < hf < 3.0 


eej and fifij 










-0.7 < hf < 0.7 


D0 Z 1 ~> [ 


142 


la 


14 


500 GeV 


-1.8 < h% < 1.8 


eej and fxp/y 










-0.5 < hf < 0.5 


D0 Zj — > 


143 


la 


13 


750 GeV 


-0.4 < fif < 0.4 


W7 










-0.06 < h% < 0.06 


D0 Z 7 -► [ 


144 


lb 


97 


500 GeV 


-1.31 < ft,f < 1.31 


ee7 and fifij 










-0.26 < hf < 0.26 


D0 Z 7 -► [ 


144 


lb 


97 


750 GeV 


-0.69 < h§ < 0.69 


ee7 and /U//7 










-0.08 < hf < 0.08 


D0 Combined [ 


144 


la + lb 




750 GeV 


-0.36 < /if < 0.36 
-0.05 < hf < 0.05 



4.4.1. W~f 

The integrated luminosity used in D0's Run la + 
Run lb W"f analysis was 6.72 times larger than the 
Run la sample alone. From that we expect the com- 
bined anomalous coupling limits to be (1/6. 72) 1 / 4 = 
0.62 as large as the Run la limits. Scaling the Run la 
results in Table |?], we would expect limits —1.0 < 
Ak 7 < 1.1 and -0.38 < A 7 < 0.38. Instead, from 
Table fj], the result was equivalent to a scaling of ~ 0.52. 
That does not seem very different, but it is, for it 
corresponds to the equivalent of a factor of two more 
luminosity. The difference is attributed to an improve- 
ment in technique, the use of a three-body transverse 
mass criteria to remove events where the photon was 
radiated from a charged final state lepton. It's hard 
to predict improvement techniques because, if such 
improvements were a priori known, they would most 
likely have been applied. 

4.4.2. WW/WZ -> ei/jj 

We expect the Run la + Run lb limits to scale by 
0.61 from the ratio of the integrated luminosities. Con- 
sulting Table [?], we find that this is essentially right on 
the nose for the Aff = 1-5 TeV limits. The Run la 
+ lb Aff = 2 TcV limits, not available because of 
unitarity constraints in the Run la sample, represent 
a slight (~ 10%) improvement over the Run la + lb 
Aff = 1.5 TeV results. 

4.4.3. WW -> dileptons 

We expect the Run la + Run lb limits to scale 
by 0.62 from the ratio of the integrated luminosities. 
However, that is not what happened. An important 



improvement in the technique, namely a 2-D fit to the 
lepton Et spectrum, plus the subsequent increase in 
the allowed form factor, allowed the combined results 
to be almost a factor of four better. Here is a case 
where we have already challenged the Run II limits 
predicted by the TeV_2000 report. 

4.4.4. Zj — > eej and fifi'y 

Again, by the fourth-root rule, we expect the Run lb 
limits to scale by 0.62. For Aff = 500 GeV, we found 
that the /if results scaled by 0.72 and the hf by 0.52, 
averaging out to a scale factor of 0.62. But, because 
of the strong dependence on the form factor scale, the 
results at Aff = 750 GeV are 3.5 and 5 times better 
for /if and hf , respectively. 

4.5. Expectations for Run II Anomalous Cou- 
pling Limits 

Having probed the usefulness and limitations of our 
scaling formula, we apply it to the Run I analyses 
to determine the limits that might be attained with 
2 fb _1 . Of course, any projections for anomalous cou- 
pling limits are merely sensitivity estimates. Improve- 
ments in technique, such as multi-dimensiona l fits , or 
using clever projection techniques (see Ref. [ 145 1 for 
an example) may yield more stringent limits. 

For the W"f and WW/WZ analyses we will use an 
integrated luminosity scale factor (2000/100) 1 / 4 = 2.1. 
The slight improvement from the w 10% increase in 
cross section available should the Tevatron operate at 
center of mass energy 2000 GeV is ignored. Increasing 
the form factor scale from 1.5 TeV to 2 TeV strength- 
ens the limits by about 10%. The WZ — > trilcptons 



39 



Table 9 

95% confidence level WWj and WW Z anomalous coupling limits that might be achieved by D0 or CDF in Run II. 



A Tl £i 1 VQl Q 


1 j l l l 1 1 . 1 UU ( 




A C Limit fQ5% CL") 

v_/ . J_J11J.11 L 1 J<J / U V / 1 J I 


W^y — » 


2000 


1.5 TeV 


-0.14 < A < 0.14 


6w~f and fJ>i // y 




Equal Couplings 


— U.44 f__ ^-iK f__ U.44 


W-y 


2000 


2.0 TcV 


-0.12 < A < 0.12 


evj and fiwy 




Equal Couplings 


-0.40 <Ak< 0.40 


WW and WZ -> 


2000 


2.0 TeV 


-0.16 < A < 0.17 






Equal Couplings 


-0.20 < An < 0.28 


VLZ — > trileptons 


2000 


2.0 TeV 


-0.2 < X z < 0.2 
-0.3 < A 5 f < 0.3 


Combined 


2000 


2.0 TeV 


-0.086 < A < 0.090 


(per experiment) 




Equal Couplings 


-0.12 < Ak < 0.19 



analysis will improve by about a factor of 6 because 
of the increased integrated luminosity and because of 
the improvement in limit-setting technique available 
by fitting the Et spectrum of the Z bosons. Table ^ 
contains the expected results. In order to put these 
bounds into perspective, we list the most recent LEP2 



0.073 < 


A 5 f 


< 0.075, 


-0.12 < 




< 0.16, 


-0.15 < 


A^ 


< 0.01. 



(95% CL) limits from a 3-parameter fit [ 146 1 , assuming 
Akz = Agf — Ak 7 tan 2 9w and Xz = A 7 : 

(17) 
(18) 
(19) 

It should be noted that form factor effects are not 
included in the bounds obtained at LEP2. Taking 
into account the form factor behavior of anomalous 
couplings weakens the limits obtained. For a dipolc 
form factor with Aff = 2 TeV, this is a 2% effect. 

For anomalous Z77 and ZZj couplings we fore- 
cast limits which are very similar to those given in 
the Tev_2000 Report. Those predictions are based on 
1 fb _1 integrated luminosity and Aff — 1500 GeV and 
are listed in Table [l0| The bounds obtained for hYo 
(h^a) almost coincide with those found for h^ {h\ ). 
For comparison, the most recent 95% CL limits on hf 
from LEP2 are: 

-0.17 <hj < 0.08 -0.26 <hf< 0.09, (20) 

-O.lK/1] <0.10 -O.lKhf < 0.16, (21) 

-0.027 <hj< 0.041 -0.29 <hf< 0.21, (22) 

-0.12 <hf< 0.20. (23) 



-0.026 <hj< 0.022 



Only one coupling at a time is varied here. Correcting 
for form factor effects, the limits for hi 3 (^2,4) weaken 
by about 5% (7%) for Aff = 1500 GeV. 

In Run II, the Tevatron will thus be able to im- 
prove the existing bounds on anomalous Z'yV cou- 
plings mostly for hif 4 . If an integrated luminosity of 



10 fb _1 can be achieved in Run II, the limits listed 
in Tables ^ and [l^ would improve by approximately a 
factor 1.5. 

A few additional comments are in order at this point: 

4.5.1. WW/WZ ei/jj 

Note that the expected Run II anomalous coupling 
limit has nearly been ruled out by the Run I combined 
analysis measurement. Nevertheless, if the Run II com- 
bined measurement is to scale based on the increase in 
the luminosity, all of the analyses must be carried out 
again. 

4.5.2. Z7 — > ee7 and fifij 

Scaling the Run I yield, totalling 29 eej+^^j candi- 
dates at D0, by the increase in luminosity, one expects 
about 600 ee7 + /i/i7 events per experiment. The QCD 
background and final state radiation background will 
be reduced compared to Run I through the application 
of a di-lepton invariant mass criteria, reducing the sam- 
ples to ~ 250 events. The Run II data will settle once 
and for all, whether there is a bump in the Zj invariant 
mass spectrum, as is not very strongly suggested by the 
Run I data. That is, unless a new one crops up. 

4.5.3. Z"f — > v 

This is a powerful tool for studying Z77 and ZZ"f 
couplings because of the large (20%) branching frac- 
tion for Z — > vv. However, because the Z boson is 
undetectable, there aren't any other kinematic han- 
dles. There are common backgrounds which produce 
the same signature as the signal: a photon recoiling 
against missing transverse energy (neutrinos). In or- 
der to reduce the backgrounds, a higher EJ, cut is 
used than in the other Zj analyses (40 GeV instead 
of 7 - 10 GeV). Understanding the normalization of 
the background from cosmic ray muons that happened 
to deposit energy in the calorimeter in such a way as to 
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Table 10 



95% confidence level Z~/V anomalous coupling limits that might be achieved by D0 or CDF experiments in Run II. 



Analysis 


Lum. (pb 1 ) 


A FF 


A.C. Limit (95% CL) 


eej 


1000 


1.5 TeV 


-0.105 < /4o < °- 105 
-0.0064 < hXo < 0.0064 




1000 


1.5 TeV 


-0.038 < hj < 0.038 
-0.0027 < h\ Q < 0.0027 



mimic a photon was the main difficulty in this analysis 
in Run I. This background should be more tractable in 
Run II using the new central and forward preshower 
detectors and a technique s imilar to that described in 
D0's Run la publications [ ^39[ [L44 |. 



W polarization at Tevatron production energies. CDF 
saw a hint [ 148 1 of the radiation zero in Run lb but 



4.6. New Directions in Di-boson Production for 
Run II 

Besides improving limits from final states analyzed 
in Run I, a number of new channels will become acces- 
sible in Run II, either due to the increased data sample, 
or because of detector improvements. In addition, it 
will be possible to search for the so-called "radiation 
zero" in Wj production. In this subsection, we first 
briefly describe the search for the radiation zero and 
the prospects for using di-boson final states involving 
6-quarks. This is followed by a somewhat more detailed 
analysis of ZZ production, the main new di-boson 
channel which will become accessible in Run II. 

4.6.1. Radiation Zero in W'y Production 

W-f production is of special interest because of 
the "radiation zero" in the helicity amplitudes 
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The Tev_2000 Report describes the situation very elo- 
quently and completely. The SM helicity amplitudes 
of the process qiq~2 — > W^j vanish for 



cosfl* 



Q1+Q2 
Qi - Q-2 



±- 



(24) 



where 9* is the scattering angle of the W boson with 
respect to the quark (q± ) direction in the W 7 rest frame 
and Qi and Q2 are the quark and anti-quark elec- 
tric charges normalized by the proton electric charge. 
Anomalous couplings destroy the radiation zero as 
do higher-order QCD corrections, backgrounds, finite 
W^-width effects, and events where the photon is radi- 
ated from the charged lepton instead of the W . 

The trick in reconstructing 9* is in determining the 
parton center-of-mass frame because there are two so- 
lutions for the z-component of the neutrino momen- 
tum. CDF discussed [ 148 ] a possible solution in se- 
lecting the minimum of the Vziy) solutions for W _ 7 
and the maximum of the Pziy) solutions for W + j. 
This is correct 73% of the time because of the high 



the signal was not definitive. 

The twofold ambiguity in reconstructing cos 9* can 
be avoided by studying rapidity correlations such as 
Ay (7,^) = 2/(7) — y{£), which manifests the radiation 
zero at Ay(j,£) w -0.3 



The Tevatron is operating at the ideal energy for 
observing the radiation zero because the zero is not 
smeared out by NLO processes expected from VF7 pro- 
duction at higher energy accelerators. Fig. [56| shows 
the Ay(-f,£) distribution, together with the statistical 
errors for 1 fb _1 . The radiation zero will be observed 
in Run II if it is there. 

4.6.2. WZ -> ivbb 

This channel has not been studied in Run I. It will 
be examined very closely in Run II because it is a 
background in the search for associated Higgs boson 
production (W + H° where W — > Iv and H° — > bb). 
The SM cross section for WZ production, including 
NLO QCD corrections, is about 3.7 pb. The branch- 
ing fraction for Z — > bb is ~ 15%. This is 2.5 times 
as much as that of Z — > /x/z and Z — > ee combined. 
So we can expect about 250 e^bb + iivbb events per 
experiment, not counting acceptance, lepton ID, and 
b-tagging efficiencies, which can be expected to amount 
to ~ 0.20 for such a final state. A cut on bb invariant 
mass will reduce the W+jets background. 

This is ripe for an anomalous coupling analysis. To 
produce WW Z anomalous coupling limits, one can fit 
the Et spectrum of the W boson and of the final-state 
lepton, such as was done in the Run I WW/WZ — > 
£i/+jets analyses. 

4.6.3. Z 7 -> bbr 

Scaling the Run I yield, totalling 29 eej + fifij can- 
didates at D0, by the increase in luminosity and a 
factor for the larger branching ratio of Z — > bb to 
Z — > ee(fifi), one might expect about 1000 Zj — > bbj 
events. Background from 7_jj and three jet events where 
a jet mimics a photon are larger than the signal and 
may constrain this to a limit-setting analysis (see also 
Sec. |3.3.5|). 
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Figure 56. Simulation of the SM Ay (7, e) distribution 
in — > W + j — > e + ^ T 7 in the D0 detector (from 



Ref. [ 77 ). In addition to the standard^, electron and 



photon pt and rapidity cuts, a AR(j, £) > 0.7 cut, and 
a cluster transverse mass cut of ?Bt(^7;pV) > 90 GeV 
are imposed (to reduce the W — > ez^y background). 
The statistical error bars for an integrated luminosity 
of 1 fb _1 are also shown. 



4.6.4. ZZ Production 

For Run I, ZZ production has not been analyzed. 
The total cross section for pp — > at -y/s = 2 TeV, in- 
cluding NLO QCD corrections is approximately 1.5 pb. 
For an integrated luminosity of 2 fb _1 one thus expects 
a few ZZ — > if £^£2^2 (£t> ^2 = e, jtt) events, if real- 
istic lepton pt and pseudo-rapidity cuts are imposed. 
Larger event rates are expected for ZZ — > £ + £~Dv, 
ZZ — > and ZZ — > Pi/jj. These channels, how- 

ever, suffer from non-trivial background contributions. 
In this subsection we briefly discuss the signals of 
anomalous ZZV couplings in the four channels, and 
derive sensitivity bounds on f\ and /g which one 
expects to achieve with 2 fb _1 (Run Ha) and 10 fb _1 
(Run lib). More details will be given elsewhere 
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The results reported here are based on a tree level 
calculation of ZZ production in the double pole ap- 
proximation. Timelike photon exchange and the de- 
cays of the Z bosons, including full decay correlations 
and finite Z width effects, are taken into account in the 
calculation. To simulate detector response, we impose 
the following transverse momentum, pseudo-rapidity 



and separation cuts: 

p T {£) > 15 GeV, 
PT (j) > 20 GeV, 
AR(£j) > 0.6, 



W)\ < 2.5, 
1*7(7)1 < 2.5, 
AR(jj) > 0.6. 



(25) 
(26) 
(27) 



In the ZZ — * i/z/jj case, Eq. ( f25|) is replaced by a 
charged lepton veto 



p T {£) < 10 GeV 



for 



\ V {£)\ < 2.5. (28) 



In addition to the cuts imposed on the leptons and jets, 
we require 



p T > 20 GeV 



< 20 GeV 
> 60 GeV 



for 
for 
for 



ZZ 
ZZ 
ZZ 



£ + ri> T , 



and 



76 GeV < m{££) < 106 GeV, 
76 GeV < m(jj) < 106 GeV. 



(29) 
(30) 
(31) 

(32) 
(33) 



Finally, in the ZZ — * £ + £~vv case, we require that 
the angle in the transverse plane between a charged 
lepton and the missing transverse momentum is be- 
tween 20° and 160° if the missing pr is j) T < 50 GeV. 
This suppresses backgrounds from bb production and 
Z — > t + t~ decays to a negligible level. 

Uncertainties in the energy measurements are taken 
into account in the numerical simulations by Gaus- 
sian smearing of the particle momenta according to 
the resolutions of the CDF II detector. For the form 
factor we use the form of Eq. m) with n = 3 and 
A FF = 750 GeV. We use the CTEQ4L parton distribu- 
tion functions with Q 2 = Af§ . Unless stated otherwise, 
only one ZZV coupling at a time is chosen different 
from its zero SM value. For simplicity, we only consider 
real ZZV couplings. 



ZZ — > 4 leptons 

Similar to the WWV and Z7U couplings, the ef- 
fects of anomalous ZZV couplings are enhanced at 
large energies. A typical signal of nonstandard ZZZ 
and ZZ7 couplings thus will be a broad increase in 
the ZZ invariant mass distribution, the Z transverse 
momentum distribution and the px distribution of the 
Z decay leptons. This is illustrated in Fig. ^7| for the 
Pt{Z) and the pt{^) distributions in pp — * ZZ — > 
e + e - /i + /i~. Results are shown for the SM, ff = 0.3 
and /j = —0.3. Terms proportional to fY and fY 
in the matrix elements have identical high energy be- 
havior. Differences in the differential cross sections at 
high energies between ZZZ and ZZ7 couplings are 
thus controlled by the Zff and 7// couplings, and 
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Figure 57. The pt(Z) and pt{p) distributions in pp 
non-standard ZZV couplings. 



ZZ 



' ' fj,~r ji in the SM and in the presence of 



by the parton distribution functions. At the Tevatron 
these result in differential cross sections which differ by 
only a few percent for a > M§ if \ ff\ = (i — 4, 5). 
Slightly larger differences are observed at intermediate 
energies and transverse momenta. Since fY violate CP 
conservation, terms in the helicity amplitudes propor- 
tional to those couplings do not interfere with the SM 
terms as long as fY is real. Cross sections thus are 
independent of the sign of fY ■ 

To distinguish fY and /g , and to determine the sign 
of fY, the AR(£+£r) and A<f>(£f£r) (i = 1,2) distri- 
butions may be helpful, if deviations from the SM pre- 
dictions should be found in the pr or the mzz differ- 
ential cross sections. Fig. |5^ shows the AR(p + p~) and 
A$(^ + distributions for pp — > ZZ — > e + e~ fi + p,~~ 
in the SM and for non-standard ZZZ couplings. The 
shape of the distributions for non-zero fY , fY > an( i 
fY < are quite different. Similar results are obtained 
for the corresponding distributions of the e + e~ pair, 
and for the ZZj couplings fj 5 . 

Anomalous couplings mostly affect the cross section 
at large if-boson transverse momentum. Due to the 
Lorentz boost, the relative opening angle between the 
leptons originating from the Z decay decreases with 
increasing px- The deviations due to non-standard 



ZZV couplings in the AR(£f£r) an d A$(£f£r) distri- 
butions thus are therefore concentrated at rather small 
values. The SM AR(£f£r) and A<f>(£f£r) differential 
cross sections are dominated by the threshold region, 
vl w 2mz, where the Z boson momenta are small 
and the decay leptons tend to be back-to-back, i.e. 
the distributions are strongly peaked at AR s» 3 and 
A$ = 180°. 

Using the AR(£f£r) and A$(£f£r) distributions, 
it may be possible to distinguish fY and fY and to 
determine the sign of fY , provided a sufficient number 
of events are observed. 



ZZ -> £+£-vv 

In contrast to the ZZ — > 4 leptons mode which is 
almost background free, there are several potentially 
important background processes if one of the two Z 
bosons decays into neutrinos. The advantage of the 
ZZ — > l + l~vv channel is its larger branching fraction. 
Summing over the three neutrino species, the number 
of ZZ — > l + l~vv signal events is about a factor 6 larger 
than the number of ZZ — > 4 leptons events. 

The most important background processes con- 
tributing to the ZZ — > l + t~vv channel are ti — » 
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W + W bb, standard electroweak W + W + X produc- 
tion with W+W- -> l+vl-v, and Z(-> + 1 jet 
production with the jet rapidity outside the range cov- 
ered by the detector and thus faking missing px- Our 
results for signal and backgrounds are summarized in 
Fig. |59| for the ZZ — > e + e~vv case. 

The two most important backgrounds are ti and 
W + W~ + X production. If no additional cuts are 
imposed to suppress the ti background, its differential 
cross section is larger than the SM signal for e+e~ 
transverse momenta as large as 200 GeV, and may 
thus reduce the sensitivity to anomalous ZZV cou- 
plings (see Fig. ^9|a). Requiring that no jets with 
Pr(j) > 20 GeV and \r)(j)\ < 3.5 are present al- 
most completely eliminates the ti background. It also 
reduces the W + W~ + X background at large trans- 
verse momenta. As shown in Fig. |55[ NLO QCD 
corrections strongly affect the px(e + e~) distribution in 
W + W~ — > e + ve~v. The enhancement at large px is 
mostly due to real emission diagrams, leading to events 
which contain a hard jet. 

To calculate the Z + 1 jet background, we have 
assumed that jets with a rapidity \rj(J)\ > 3.5 are 
misidentified as j> T . With this rather conservative as- 
sumption, the Z+l jet background is much larger than 
the ZZ signal at small transverse momenta. Due to 
kinematical constraints, however, it drops rapidly with 
Px- Since non-standard ZZV couplings lead to large 
deviations from the SM only at high transverse mo- 
mentum, essentially no sensitivity is lost by requiring 
p T (£ + £~) > 40 GeV when testing for ZZV couplings. 

ZZ -> £+£~jj and ZZ -> vvjj 

The ZZ — > ^ + ^~jj and ZZ — > vvjj channels have 
larger branching ratios than the ZZ — > 4 leptons and 
the ZZ — ► £ + t~vv channels, but also much higher 
backgrounds. The main background sources are QCD 
Z + 2 jet production and W+Z production with the 
W decaying into two jets. The px{£ + £~) distribution 
for ZZ — > £ + £~]j is shown in Fig. |6C|a. Fig. |60| b shows 
the px(jj) distribution for ZZ — > i>vjj. In each case 
we display the SM cross section together with the two 
main backgrounds, Zjj and W+Z production. We also 
show the ZZ cross section for ff = 0.3. 

The p'rp > 60 GeV cut imposed in the ZZ — > ^ T jj 
case helps to suppress the bb and Z — ► t + t - back- 
grounds. The "kink" in the WZ and ZZ differen- 
tial cross sections at px ~ 250 GeV is due to the 
AR(jj) > 0.6 cut which becomes effective only at 
sufficiently high transverse momenta. The W+Z dif- 
ferential cross section is very similar to that of the 
SM signal over most of the px range considered. The 
j> T < 20 GeV cut imposed in the ZZ — > £ + £~]j channel 



effectively eliminates the ti — > £ + v£~i>jj background. 
The charged lepton veto (see Eq. (p^)) required in the 
ZZ — > ^ r jj case rejects backgrounds from ti produc- 
tion, W — + lu, and Z — * decays. The Zjj back- 
ground is uniformly about a factor 10 larger than the 
SM ZZ signal. It will therefore be very difficult to ob- 
serve ZZ production in the semi-hadronic channels, if 
the SM prediction is correct. However, for sufficiently 
large anomalous ZZV couplings, the ZZ cross section 
exceeds the background at large transverse momenta. 
ZZ — > £ + £~]] and ZZ — > vv]] therefore may still be 
useful in obtaining limits on the ZZV couplings, sim- 
ilar to the semi-hadronic WW and WZ channels used 
by CDF and D0 in Run I to extract limits on the 
WWV couplings. 

Sensitivity Bounds 

In order to derive sensitivity limits for anomalous 
ZZV couplings which one can hope to achieve in 
Run II, we use the px{£ + £~) distribution for ZZ — » 
4 leptons, ZZ -> £+£rDv and ZZ -> ^ - jj. For 
the ZZ — > Dv']] channel we use the px{]\) distribu- 
tion. Other distributions, such as the ZZ invariant 
mass distribution (useful only for ZZ — > 4 leptons), 
or the maximum or minimum transverse momenta of 
the charged leptons or jets, yield similar results. In 
deriving our sensitivity limits, we combine channels 
with electrons and muons in the final state. 

We calculate 95% confidence level (CL) limits per- 
forming a x 2 test. The statistical significance is calcu- 
lated by splitting the px distribution into a number of 
bins, each with more than five events typically. In each 
bin the Poisson statistics is approximated by a Gaus- 
sian distribution. In order to derive realistic limits, 
we allow for a normalization uncertainty of 30% of the 
SM cross section. Backgrounds in the ZZ — > £ + £~vv, 
ZZ — > ^ + ^ _ jj and ZZ — > vv]] channels are included 
in our calculation. In the ZZ — * £ + £~vv case we 
assume that a jet veto has been imposed to reduce 
the ti background and require px(£ + £~) > 40 GeV to 
eliminate the Z + 1 jet background. As before, we use 
a form factor of the form of Eq. (H) with n = 3 and 
App = 750 GeV. Non-negligible interference effects 
are found between ff and fj : and between /g and 
f2 ■ As a result, different anomalous contributions to 
the hclicity amplitudes may cancel partially, resulting 
in weaker bounds than if only one coupling at a time 
is allowed to deviate from its SM value. 

In Table [H] we display sensitivity limits for the 
Tevatron and integrated luminosities of 2 fb _1 and 
10 fb _1 , taking into account the correlations between 
ff and fj , and between f§ and /g. No limits for the 
ZZ — > 4 leptons case with 2 fb _1 are given. The lim- 
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Figure 60. Transverse momentum distribution of a) the £ + £ pair in pp — > ZZ — > £ + £ jj, and b) of the jet pair in 
ZZ — > Dvjj at the Tevatron, together with the differential cross sections from ti, W^Z and Z + 2 jet production. 



ited number of events in this case does not allow for an 
analysis of the px(£ + £~) distribution using the method 
chosen here. The bounds obtained from ZZ — » £ + £~vv 
and ZZ — > vv]] are quite similar. The cross section 
for ZZ — > Pi/jj is about a factor 10 larger than that 
for ZZ — > £ + £~vv, however, the large background 
from Zjj production considerably limits the sensitiv- 
ity to ZZV couplings for ZZ — > vv]]. The limits 
from the ZZ — > ^ + ^~jj and — > 4 leptons channels 
are about a factor 1.5 and 2 weaker than those from 
ZZ — > £ + £~Du and ZZ — > wjj. 

We have not made any attempt to combine the limits 
from different channels. From Table [ll] it is clear that 
this would result in a significant improvement of the 
bounds. 

The sensitivity limits which can be achieved at the 
Tevatron in Run II should be compared with the 
bounds from recent measurements at LEP2 [ |14£ 



\.fi 



< 0.49 \fj\ < 0.82 (34) 

\f 5 z \ < 1.1 |#| < 1.1 (35) 

Only one coupling at a time is varied here. (The LEP2 
limits do not contain any form factor effects. For the 
form and scale chosen here, form factor effects weaken 
the limits by about 20%.) In Run II, CDF and D0 will 
be able to improve these bounds by at least a factor 4 



to 8. 

4.6.5. Measuring Form Factors 

The limits on anomalous WWV ', ZjV and ZZV 
couplings all depend on the power, n, and the scale, 
Aff, °f the form factor. These parameters are a pri- 
ori unknown. In Ref. [ |117| it was pointed out that 
in final states without missing transverse momentum 
one can in principle determine the form factor by mea- 
suring the y/i distribution. For W^y production, the 
longitudinal momentum of the neutrino can only be 
reconstructed with a twofold ambiguity. Selecting the 
minimum of the two reconstructed values of s, a si mila r 
measurement can be performed in the Wj case 
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Alternativel y, th e photon px distribution can be used. 

In Ref. [ 127 1 a detailed study of the method was 
performed for W(— > ev, pis)"/ production at the LHC. 
Assuming Aq = 0.025, n = 2, A F f = 2 TeV and an 
integrated luminosity of 300 fb _1 , AJJ and Aff were 
reconstructed using a binned maximum likelihood fit 
to the Pt(i) distribution. The power of the form 
factor, n = 2, was not varied and no detector simu- 
lation was included in the study. The reconstructed 
coupling and form factor scale were determined to 
\2 = 0.0295 ± 0.0022 and A FF = 1.67 ± 0.22 TeV, 
ie. they can be measured with a relative precision 
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Tabic 11 

Sensitivities achievable at 95% CL for the anomalous ZZV couplings in pp — > ZZ — ► 4 leptons, pp — > ZZ — » l+l~vv, 
pp — ► — > and — ► ZZ — ► Pz/jj at the Tevatron a) for an integrated luminosity of 2 fb _1 , and b) for an 

integrated luminosity of 10 fb _1 . 



a) JCdt = 2 fbr 1 

coupling ZZ -> 4 leptons ZZ -> £ + ^"z7z/ ZZ -> ZZ -> Pi/jj 



f z 

J40 


- 


+0.169 
-0.169 


+0.219 
-0.218 


+0.156 
-0.155 


flo 


- 


+0.175 
-0.174 


+0.222 
-0.221 


+0.157 
-0.157 


fib 


- 


+0.171 
-0.204 


+0.220 
-0.244 


+0.157 
-0.179 


fhO 


- 


+0.184 
-0.202 


+0.229 
-0.241 


+0.166 
-0.174 


coupling 


ZZ — > 4 leptons 


b) JCdt = 10 fb- 1 
ZZ -» £+^"z7z/ 


ZZ -> 


ZZ - Pi/jj 


fio 


+0.180 
-0.179 


+0.097 
-0.097 


+0.146 
-0.145 


+0.104 
-0.103 


flo 


+0.185 
-0.185 


+0.100 
-0.099 


+0.148 
-0.147 


+0.104 
-0.104 


f z 

J 50 


+0.178 
-0.216 


+0.092 
-0.120 


+0.144 
-0.167 


+0.102 
-0.124 


fbO 


+0.192 
-0.213 


+0.103 
-0.115 


+0.151 
-0.163 


+0.109 
-0.118 



of about 10 - 15%. The central values of the recon- 
structed parameters differ by about 20% from the input 
parameters. Including detector response, and treating 
the form factor power n as an additional free parameter 
is expected to increase the relative error. 

The study demonstrates that, due to the restricted 
number of events in each bin, the method will not pro- 
duce competitive limits. However, if non-zero anoma- 
lous couplings are observed, the method may be use- 
ful in determining the shape of the form factor which 
provides indirect information on the dynamics of the 
underlying new physics. 

At the Tevatron, the limits on the WWV couplings 
exhibit only a modest dependence on n and Aff- Di- 
rect measurement of the form factor will thus be very 
difficult for these couplings. The situation is more 
promising for Z"fV and ZZV couplings, where the 
sensitivity bounds depend more strongly on the form 
factor parameters. 

Acknowledgements 

We acknowledge contributions from the other mem- 
bers of the working group including P. Aurcnchc, 



C. Balazs, L. Dixon, S. Eno, G. Gomez, N. Kidon- 
akis, P. Nason, J. Qiu, and G. Sterman. We also 
thank T. Ferbel, M. Fontannaz, G. Ginthcr, B. Kniehl, 
A. D. Martin, R. Roberts, P. Slattery, J. Stirling, 
R. Thorne and W.K. Tung. We also wish to thank 
A. Beretvas, C. Bigongiari, J. Campbell, O. Lobban, 
and W. Wester for their help. Special thanks go to 
M. Mangano for some enlightening discussions. 

REFERENCES 

1. R. Brock et al, Report of the Precision Measure- 
ment Working Group. 

2. L. dcBarbaro et al, Report of the Working Group 
on Parton Distribution Functions. 

3. F. Abe et al. (CDF Collaboration), Phys. Rev. 
D60, 092003 (1999). 

4. B. Bailey, E L. Berger, and L. E. Gordon, 
Phys.Rev. D54, 1896 (1996). 

5. P. Aurcnchc et al, Phys. Lett. B140, 87 (1984). 

6. E.L. Berger, E. Braaten, and R.D. Field, Nucl. 
Phys. B239, 52 (1984). 

7. E.L. Berger and J.W. Qiu, Phys. Rev. D44,2002 
(1991). 



47 



8. B. Bailey, J. F. Owens, and J. Ohncmus, Phys. 
Rev. D46, 2018 (1992). 

9. T. Binoth et al, Eur. Phys. J. C16, 311 (2000). 

10. C. Balazs, E.L. Berger, S. Mrenna, and C.-P. 
Yuan, Phys. Rev. D57, 6934 (1998). 

11. P. Chiappetta, R. Fergani, and J. Guillet, Phys. 
Lett. B348, 646 (1995). 

12. S. Catani et al, JHEP 9903, 025 (1999). 

13. N. Kidonakis and J. F. Owens, Phys. Rev. D61, 
094004 (2000). 

14. E. Laenen, G. Oderda, and G. Sterman, Phys. 
Lett. B438, 173 (1998). 

15. H.-L. Lai and H. nan Li, Phys. Rev. D58, 114020 

(1998) . 

16. H. nan Li, Phys. Lett. B454, 328 (1999). 

17. E. Laenen, G. Sterman, and W. Vogelsang, Phys. 
Rev. Lett. 84, 4296 (2000). 

18. J. Huston et al. (CTEQ Collaboration), Phys. 
Rev. D58, 114034 (1998). 

19. H. L. Lai et al. (CTEQ Collaboration), Phys. Rev. 
D55, 1280 (1997). 

20. M. Chick, E. Reya, and A. Vogt, Z. Phys. C53, 
127 (1992). 

21. A. D. Martin et al, Eur. Phys. J. C4, 463 (1998). 

22. H. Baer and M. Reno, Phys. Rev. D54, 2017 
(1996). 

23. J. Huston et al. (CTEQ Collaboration), Phys. 
Rev. D51, 6139 (1995). 

24. L. Apanasevich et al, Phys. Rev. D59, 074007 

(1999) . 

25. P. Aurenche et al, Eur. Phys. J. C9, 107 (1999). 

26. P. Aurenche et al, Eur. Phys. J. C13, 347 (2000). 

27. M.A. Kimber, A.D. Martin, and M.G. Ryskin, 
Eur. Phys. J. C12, 655 (2000). 

28. C. Balazs, J. Huston, and I. Puljak, hep 



ph/0002032| (February 2000), submitted to Phys. 



Rev. D. 

29. C. Balazs, Ph.D. thesis, Michigan State Univer- 
sity, 1999. 

30. S. Catani et al, |hep-ph/0005114| (May 2000), to 
appear in the Proceedings of the Workshop on 
Physics at TcV Colliders (Les Houches, 1999). 

31. M. Begel, Nucl. Phys. B (Proc. Supp.) 79, 244 
(1999); Ph.D. thesis, University of Rochester, 
1999. 

32. E. Bonvin et al. (WA70 Collaboration), Phys. 
Lett. B236, 523 (1990). 

33. F. Abe et al (CDF Collaboration), Phys. Rev. 
Lett. 70, 2232 (1993). 

34. R. Barate et al. (WAll Collaboration), Phys. Rev. 
Lett. 21, 1541 (1979); M. J. Corden et al (WA12 
Collaboration), Phys. Lett. B76, 226 (1978); 
A. L. S. Angelis et al. (R108 Collaboration), Phys. 
Lett. B97, 163 (1980); D. Antreasyan et al. (R209 



Collaboration), Phys. Rev. Lett. 47, 12 (1981); 
J. Badier et al (NA3 Collaboration), Phys. Lett. 
B117, 372 (1982); A. S. Ito et al (E288 Collabora- 
tion), Phys Rev. D23, 604 (1981); D. Antreasyan 
et al (E300 Collaboration), Phys. Rev. Lett. 39, 
906 (1977); M. D. Corcoran et al (E395 Collabo- 
ration), Phys. Rev. D21, 641 (1980); K. J. Ander- 
son et al. (E444 Collaboration), Phys. Rev. Lett. 
42, 944 (1979); E. Anassontzis et al (E537 Collab- 
oration), Phys. Rev. D38, 1377 (1988); G. Moreno 
et al (E605 Collaboration), Phys. Rev. D43, 2815 
(1991); S. Palestini et al. (E615 Collaboration), 
Phys. Rev. Lett. 55, 2649 (1985); J. S. Conway 
et al, Phys. Rev. D39, 92 (1989); J. Alspector 
et al. (E687 Collaboration), Phys. Lett. B81, 397 
(1979). 

35. A. L. S. Angelis et al (R108 Collaboration), Phys. 
Lett. B97, 163 (1980); A. G. Clark et al (R702 
Collaboration), Nucl. Phys. B160, 397 (1979). 

36. W. Chen, Ph.D. thesis, State University of New 
York at Stony Brook, 1997. 

37. H.-U. Bengtsson and T. Sjostrand, Comput. Phys. 
Commun. 46, 43 (1987), [Pythia v5.7]. 

38. J. F. Owens, Rev. Mod. Phys. 59, 465 (1987). 

39. A. P. Contogouris, S. Papadopoulos, and C. Pa- 
pavassiliou, Nucl. Phys. B179, 461 (1981). 

40. A. P. Contogouris et al, Phys. Rev. D32, 1134 
(1985). 

41. M. Fontannaz and D. Schiff, Nucl. Phys. B132, 
457 (1978). 

42. R. P. Feynman, R. D. Field, and G. C. Fox, Phys. 
Rev. D18, 3320 (1978). 

43. J. Collins and D. Soper, Nucl. Phys. B193, 381 
(1981); B213, 545(E) (1983); B197, 446 (1982); 
J. Collins, D. Soper, G. Sterman, Phys. Lett. 
B109, 388 (1982). 

44. J.C. Collins, D.E. Soper, and G. Sterman, Nucl. 
Phys. B250, 199 (1985). 

45. L. Apanasevich et al (E706 Collaboration), Phys. 
Rev. Lett. 81, 2642 (1998). 

46. M. Bonesini et al. (WA70 Collaboration), Z. Phys. 
C37, 535 (1988). 

47. M. Bonesini et al. (WA70 Collaboration), Z. Phys. 
C38, 371 (1988). 

48. G. Ballocchi et al (UA6 Collaboration), Phys. 
Lett. B436, 222 (1998). 

49. E. Bonvin et al (WA70 Collaboration), Z. Phys. 
C41, 591 (1989). 

50. L. Apanasevich, Ph.D. thesis, Michigan State Uni- 
versity, 2000. 

51. M. Bonesini et al. (WA70 Collaboration), Z. Phys. 
C37, 39 (1987). 

52. J. Antille et al. (UA6 Collaboration), Phys. Lett. 
B194, 568 (1987). 



48 



Collaboration) , 
Collaboration), 
Phys. 
79, 



53. E. Anassontzis et al. (R806 
Z. Phys. C13, 277 (1982). 

54. C. Kourkoumelis et al. (R806 
Z. Phys. C5, 95 (1980). 

55. M. Diakonou et al. (R806 Collaboration) 
Lett. B91, 296 (1980). 

56. S. Kuhlmann, Nucl. Phys. B (Proc. Suppl 
241 (1999). 

57. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
Lett. 84, 2786 (2000). 

58. R. Ansari et al. (UA2 Collaboration), Z. Phys. 
C41, 395 (1988). 

59. J. Alitti et al. (UA2 Collaboration), Phys. Lett. 
B288, 386 (1992). 

60. S. Abachi et al. (D0 Collaboration), FERMILAB- 
PUB-96-357-E. 

61. V. Zutshi, D0 note in preparation. 

62. Y. Gershtein, D0 note in preparation. 

63. H. L. Lai et al. (CTEQ Collaboration), Phys. Rev. 
D51, 4763 (1995). 

64. H. L. Lai et al. (CTEQ Collaboration), Eur. Phys. 
J. C12, 375 (2000). 

65. G. Sterman and W. Vogelsang, [hcp-ph/00021~32| 
(February 2000). 

66. S. Catani, M.L. Mangano, and P. Nason, hep 



ph/9806484, JHEP 9807, 024 (1998). 



67. 



68. 



69. 



70. 



71. 



72. 



73. 



74. 



75. 



76. 



77. 



78. 



79. 



hcp-ph/9706526, Nucl 



Yu.L. Dokshitzer, D.I. Dyakanov, and S.I. Troyan, 
Phys. Rep. 58, 269 (1980). 
R.K. Ellis and S. Veseli, ~ 
Phys. B511, 649 (1998). 

W. Vogelsang and A. Vogt, |hcp-ph/9505404| , Nucl. 
Phys. B453, 334 (1995). 

H. Baer, J. Ohncmus, and J.F. Owens, Phys. Rev. 
D42, 61 (1990). 

F. Abe et al. (CDF Collaboration), Phys. Rev. 
Lett. 77, 439 (1996); a paper on the Run lb anal- 
ysis is in preparation. 

B. Abbott et al. (D0 Collaboration), Phys. Rev. 
Lett. 82, 2451 (1999). 

F. Abe et al. (CDF Collaboration), Phys. Rev. 
D57, 1359 (1998). 

E. L. Berger, Xiao-feng Guo, and Jian-wei Qiu, 



hep-ph/9512281, Phys.Rev.Lett. 76, 2234 (1996); 



hep-ph/9605324| , Ph ys. Rev. D54, 54 70 (1996) 
hep-ph/9606287], Phys 



Rev. 



P. Aurenche et al.. 

D55 1124 (1 997). 

S. Frixionc, |hcp-ph/9801442| , Phys. Lett. B429, 
369 (1998). 

D. Amidei et al., Future Electro-Weak Physics at 
the Fermilab Tevatron, edited by D. Amidei and 
R. Brock, FERMILAB-Pub-96/082 (1996). 
B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D61, 072001 (2000). 

D. Cronin-Hennessy, A. Beretvas, P.F. Derwent, 



Nucl. Instrum. Mcth. A443, 37 (2000). 

80. F. Abe et al. (CDF Collaboration), Phys. Rev. 
D50, 5550 (1994). 

81. C. Avila et al. (E811 Collaboration), Phys. Lett. 
B445, 419 (1999). 

82. N. Amos et al. (E710 Collaboration), Phys. Rev. 
Lett. 68, 2433 (1992). 

83. M. Albrow, A. Beretvas, L. Nodulman, and 
P. Giromini, CDF-PUB-4844 (1999). 

84. J. Bantly, J. Krane, D. Owen, R. Partridge, and 
L. Patcrno, Fcrmilab-TM-1995 (1997). 

85. R. Hamberg, W.L. van Neerven, and T. Matsuura, 
Nucl. Phys. B359, 343 (1991); W.L. van Neerven 
and E.B. Zijlstra, Nucl. Phys. B382, 11 (1992). 

86. For a review of resummation in transverse mo- 
mentum space, see R.K. Ellis and S. Veseli, Nucl. 
Phys. B511, 649 (1998). 

87. For a review of resummation in the Fourier con- 
jugate impact parameter space, see C. Balazs and 
C.-P. Yuan, Phys. Rev. D56, 5558 (1997). 

88. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D58, 092003 (1998). 

89. F. Abe et al. (CDF Collaboration), Phys. Rev. 
D52, 4784 (1995). 

90. C. Balazs and C.-P. Yuan, Phys. Rev. D59, 
114007 (1999). 

91. Q. Fan and A. Bodek, Proceedings, 6 th Interna- 
tional Conference on Calorimetry in High Energy 
Physics (ICCHEP 96), Rome, Italy, 8-14 June 
1996, pp. 553 - 560. 

92. T. Sjostrand, Comput. Phys. Commun. 82, 74 
(1994). 

93. E. Barberio and Z. Was, Comput. Phys. Commun. 
79, 291 (1994); E. Barberio, B. van Eijk, and Z. 
Was, ibid. 66, 115 (1991). 

94. T. Affolder et aZ.(CDF Collaboration), Phys. Rev. 
Lett. 84, 845 (2000). 

95. G.A. Ladinsky and C.-P. Yuan, Phys. Rev. D50, 
4239 (1994). 

96. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D61, 032004 (2000). 

97. U. Baur, S. Keller, and W. K. Sakumoto, Phys. 
Rev. D57, 199 (1998). 

98. F. Abe et al. (CDF Collaboration), Phys. Rev. 
D59, 052002 (1999). 

99. A. Bodek et al. for the CDF Collaboration, Con- 
straints on PDF's from W and Z Rapidity Distri- 
butions at CDF, Fermilab-Conf-99/160-E, to be 
published in Proceedings, 7 th International Con- 
ference on Deep Inelastic Scattering and QCD 
(DIS 99), DESY, Zeuthen, Germany, 19 - 24 April 
1999. 

100. G. Arnison et al. (UA1 Collaboration), Phys. Lett. 
B122, 103 (1983). 



49 



101. G. Arnison et al. (UA1 Collaboration), Phys. Lett. 
B126, 398 (1983). 

102. G. Banner et al. (UA1 Collaboration), Phys. Lett. 
B122, 476 (1983). 

103. P. Bagnaia et al. (UA1 Collaboration), Phys. Lett. 
B129, 130 (1983). 

104. J. Alitti et al. (UA2 Collaboration), Phys. Lett. 
B186, 452 (1987). 

105. J. Alitti et al. (UA2 Collaboration), Z. Phys. C49, 
17 (1991). 

106. M. Carena et al., Report of the Higgs Working 
Group, Fermilab Supersymmetry/Higgs Run II 
Workshop (1998). 

107. U. Baur and E.L. Berger, Phys. Rev. D47, 4889 
(1993). 

108. E. Mirkes, Nucl. Phys. B387, 3 (1992). 

109. J. C. Collins and D.E. Soper, Phys. Rev. D16, 
2219 (1977). 

110. M.I. Martin, Measurement of the Angular Distri- 
bution of the Electron from W — > ev decay in pp at 
yfs= 1.8 TeVas Function ofp%, PhD thesis, Uni- 
versidad de Zaragoza, Spain, 1994 (unpublished). 

111. G. Steinbruck, Measurement of the Angular Dis- 
tribution of Electrons from W Boson Decays at 
D0 , PhD thesis, University of Oklahoma, Nor- 
man, Oklahoma, 1999 (unpublished). 

112. K. Hagiwara, R. D. Peccei, D. Zeppenfeld and 
K. Hikasa, Nucl. Phys. B282, 253 (1987). 

113. U. Baur and D. Zeppenfeld, Nucl. Phys. B308, 
127 (1988). 

114. G. Gounaris et al., in Physics at LEP2, edited by 
G. Altarelli, T. Sjostrand and F. Zwirner, CERN 
96-01, 525 (1996). 

115. K. Hagiwara, S. Ishihara, R. Szalapski and D. Zep- 
penfeld, Phys. Rev. D48, 2182 (1993). 

116. U. Baur and D. Rainwater, FERMILAB-Pub- 
00/153-T, to appear in Phys. Rev. D. 

117. G.J. Gounaris, J.Layssac and F.M. Renard, Phys. 
Rev. D61, 073013 (2000). 

118. B. Mele, P. Nason and G. Ridolfi, Nucl. Phys. 
B357, 409 (1991); S. Frixione, P. Nason and 
G. Ridolfi, Nucl. Phys. B383, 3 (1992); S. Frix- 
ione, Nucl. Phys. B410, 280 (1993). 

119. U. Baur, T. Han, and J. Ohnemus, Phys. Rev. 
D48, 5140 (1993). 

120. U. Baur, T. Han, and J. Ohnemus, Phys. Rev. 
D57, 2823 (1998). 

121. U. Baur, T. Han, and J. Ohnemus, Phys. Rev. 
D53, 1098 (1996). 

122. U. Baur, T. Han, and J. Ohnemus, Phys. Rev. 
D51, 3381 (1995). 

123. J.M. Campbell and R.K. Ellis, Phys. Rev. D60, 
113006 (1999). 

124. L. Dixon, Z. Kunszt and A. Signer, Phys. Rev. 



D60, 114037 (1999). 

125. D. DeFlorian and A. Signer, Eur. Phys. J. C16, 
105 (2000). 

126. H.L. Lai et al. (CTEQ Collaboration), Phys. Rev. 
D55, 1280 (1997). 

127. S. Haywood et al., Report of the Electroweak 



Physics Working Group, hep-ph/0003275 (March 
2000), to appear in the CERN Yellow- Report of 
the 1999 CERN Workshop on Standard Model 
Physics (and more) at the LUC. 

128. M. Dobbs, private communication. 

129. F. Abe et al. (CDF Collaboration), Phys. Rev. 
Lett. 74, 1936 (1995). 

130. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 75, 1034 (1995). 

131. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 78, 3634 (1997). 

132. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 75, 1023 (1995). 

133. F. Abe et al. (CDF Collaboration), Phys. Rev. 
Lett. 78, 4537 (1997). 

134. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D58, 051101 (1998). 

135. F. Abe et al. (CDF Collaboration), Phys. Rev. 
Lett. 75, 1017 (1995). 

136. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 77, 3303 (1996). 

137. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
Lett. 79, 1441 (1997). 

138. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D60, 072002 (1999). 

139. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
D56, 6742 (1997). 

140. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D58, 031102 (1998). 

141. F. Abe et al. (CDF Collaboration), Phys. Rev. 
Lett. 74, 1941 (1995). 

142. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 75, 1028 (1995). 

143. B. Abbott et al. (D0 Collaboration), Phys. Rev. 
D57, 3817 (1998). 

144. S. Abachi et al. (D0 Collaboration), Phys. Rev. 
Lett. 78, 3640 (1997). 

145. J. Papavassiliou and K. Philippidcs, Phys. Rev. 
D60, 113007 (1999). 

146. G. Bella et al. (The LEP TGC Working Group), 
LEPEWWG/TGC/2000-01 (March 2000). 

147. RW. Brown and K.O. Mikaelian, Phys. Rev. 
D19, 922 (1979). 

148. D. Benjamin, Proceedings of the 10th Topical 
Workshop on Proton- Antiproton Collider Physics, 
Batavia, IL, May 1995, edited by R. Raja and 
J. Yoh, AIP Press (1996), p. 370. 

149. U. Baur, S. Errede and G. Landsberg, Phys. Rev. 



50 



D50, 1917 (1994). 



